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АННОТАЦИЯ. На примере г. Ижевска показано влияние извлечения сортируемых 
фракций на теплотехнические характеристики твердых бытовых отходов (ТБО). 
Оценено содержание в ГБО влаги, горючих веществ и золы в разные сезоны года, а 
также их морфологический состав. Приведены результаты расчетов теплоты сгора- 
ния ТБО. В результате сортировки происходит изменение свойств получаемого суб- 
страта ТБО, что может сказаться при их утилизации различными способами. 


ВВЕДЕНИЕ 


Неизбежным следствием развития современной цивилизации является все возрастаю- 
щее количество отходов. Лишь небольшая часть их подвергается утилизации и обез- 
вреживанию. Большая часть, попадая в окружающую среду, является источником по- 
ступления в эту среду различных токсичных веществ. 

Для выбора метода обезвреживания ТБО необходимо тщательно исследовать их 
морфологический и фракционный состав, плотность. Состав ТБО, а следовательно их 
свойства, изменяются в течение года. Летние бытовые отходы содержат значительно 
больше растительных составляющих с большим количеством влаги, в зимних отходах 
наибольшая доля золы. 

Наблюдается определенная закономерность в перспективном изменении этого 
состава - предвидится значительное увеличение количества бумаги и пластмасс, а сле- 
довательно снижение плотности и влажности, увеличение теплоты сгорания отходов. 

Основной способ обезвреживания ТБО большинства стран - складирование на 
полигонах [2, 3]. От российского города с миллионным населением в настоящее время 
поступает не городские свалки в течение года примерно 260 тыс. тонн ТБО. При сред- 
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ней высоте складирования мусора 10 м и плотности ТБО 0,2 т/м 3 ежегодно требуемая 
территория под свалки для такого города составит 13,5 га [1]. 

На свалках (полигонах), по существу, происходит захоронение ценных ресурсов, 
которые могли бы быть полезно использованы. Проблема истощения ресурсов из-за 
расточительного использования металлов и других элементов частично может быть 
решена благодаря усовершенствованию методов повторного использования и цикличе- 
ской переработки. Многие компоненты ТБО возможно переработать в полезные про- 
дукты, из отходов может быть выделено следующее вторичное сырье: бумага (картон, 
газеты), текстиль, стекло, авторезина, лом цветных и черных металлов, пластик, акку- 
муляторы и бытовые батарейки. 

Вторичной переработке (в т.ч. переработке в компост) подвергается незначи- 
тельное количество ТБО. Так по данным [2] вторичной переработке в США подверга- 
ется 13,1% всех ТБО, в ФРГ - около 15% всех ТБО, а в РФ - лишь 1,3% [2,3]. 

Для эффективного соора, переработки и удаления ГБО необходимо дальновид- 
ное региональное планирование. При составлении муниципальных программ [4,6] по 
решению проблемы ТБО необходимо использовать комплексный подход [2], который 
бы учитывал имеющиеся и перспективные возможности в регионе по переработке вто- 
ричных ресурсов. На основании этого может быть заложен тот или иной метод решения 
проблемы переработки ТБО. 

Одним из способов решения проблемы ТБО для городов с численностью насе- 
ления до 1 млн. жителей является создание линии сортировки мусора (ЛСМ). ЛСМ мо- 
гут применяться самостоятельно, либо в комплексе с другими способами обработки 
ТБО. В данной статье проведен анализ возможности использования ЛСМ на примере г. 
Ижевска. 

ЗАВИСИМОСТЬ СВОЙСТВ МУСОРА ОТ СОРТИРОВКИ НА ЛСМ 

В результате сортировки происходит изменение свойств получаемого субстрата ТБО. 
что может сказаться при их утилизации различными способами, а также при захороне- 
нии на специализированном полигоне. 

Для изучения зависимости свойств ТБО от степени сортировки использованы 
данные исследований, проведенных Академией коммунального хозяйства (АКХ им. 
К.Д.Памфилова), г. Москва [5]. Результаты этих исследований были взяты в качестве 
исходных данных и проведена оценка свойств ТБО г. Ижевска при различной степени 
извлечения отдельных фракций (в % масс), которые могут использоваться в качестве 
вторичного сырья. 

В таблице 1 приведены результаты расчета средних значений морфологического 
состава ТБО г. Ижевска при различной степени извлечения сортируемых фракций. 
Приведенная в табл.1 и далее степень извлечения фракций в 100% не имеет практиче- 
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Таблица 1 . Средний расчетный морфологический состав ТБО г. Ижевска в летний, 

осенний и зимний периоды (% масс.), в зависимости от степени извлечения 
сортируемых фракций 


Фракция 

Степень извлечения сортируемых фракций, 

% от исходного содержания 



0 

20 

40 

60 

80 

100 

Бумага 

лето 

20,4 

18,1 

15,2 

11,6 

6,7 

0 

осень 

17,6 

15,3 

12,5 

9,2 

5,1 

0 


зима 

20,0 

17,4 

14,3 

10,6 

6,0 

0 

Пище- 

лето 

35,6 

39,6 

44,4 

50,4 

58,4 

69,8 

вые 

осень 

41,1 

44,5 

48,7 

53,8 

59,8 

67,6 

отходы 

зима 

44,8 

48,9 

53,6 

59,7 

67,0 

76,2 

Садовые 

лето 

2,9 

3,2 

3,6 

4,1 

4,8 

5,7 

отходы 

осень 

8,1 

8,8 

9,6 

10,6 

11.8 

13,3 


зима 

0 

0 

0 

0 

0 

0 

Дерево 

лето 

0,9 

1,0 

1,1 

1,3 

1,5 

1,8 

осень 

2,2 

2,4 

2,6 

2,9 

3,2 

3,6 


зима 

1,7 

1,8 

2,0 

2,3 

2,5 

2,9 

Черный 

лето 

4,3 

3,8 

3,2 

2,4 

1,5 

0 

металло- 

осень 

3,6 

3,1 

2,6 

1,8 

1,0 

0 

лом 

зима 

3,0 

2,6 

2,2 

1,6 

0,9 

0 

Цветной 

лето 

1,6 

1,4 

1,2 

0,8 

0,5 

0 

металло- 

осень 

0,7 

0,6 

0,5 

0,4 

0,2 

0 

лом 

зима 

0,6 

0,5 

0,4 

0,3 

0,2 

0 

Кости 

лето 

1,9 

2,1 

2,4 

2,7 

3,1 

3,7 

живот- 

осень 

1 ,0 

М 

1,2 

1,3 

1,5 

1,6 

ных 

зима 

2,5 

2,7 

3,0 

3,3 

3,7 

4,3 

Резина, 

лето 

1,8 

1,6 

М 

1,0 

0,7 

0 

кожа 

осень 

1,4 

1,2 

0,9 

0,7 

0,4 

0 


зима 

3,0 

2,6 

2,2 

1,6 

0,9 

0 

Текстиль 

лето 

8,7 

7,8 

6,5 

5,0 

2,8 

"0 


осень 

7,0 

6,1 

5,0 

3,7 

2,0 

0 


зима 

5,5 

4,8 

3,9 

2,9 

1,6 

0 


лето 

1- 

6,8 

6,0 

5,1 

3,8 

2,3 

0 

Стекло 

осень 

2,8 

2,4 

2,0 

1,4 

0,9 

0 


зима 

5,6 

4,9 

4,1 

2,9 

1,6 

0 
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Строи- 

тельный 

мусор, 

камни 

лето 

1,4 

ГТб“ 

1,7 

2,0 

2,3 

2,7 

осень 

0,8 

0,9 

0,9 

1,0 

" 1,2 

1,3 

зима 

1,6 

1,7 

1,9 

2,1 

2,4 

2,7 

Пласт- 

масса 

лето 

5,3 

4,6 

3,9 

3,1 

1.8 

0 

осень 

6,1 

5,3 

4,4 

3,1 

1,8 

0 

зима 

3,6 

3,2 

2,6 

1,8 

1,0 

0 

Отсев 
фракций, 
менее 16 

мм 

лето 

8,3 

Г^2~ 

10,3 

11,8 

13,6 

16,3 

осень 

7,7 

8,3 

9,1 

10,1 

11,2 

12,6 

зима 

8,2 

8,9 

9,8 

10,9 

12,2 

13,9 


ского значения и показывает лишь возможное предельное значение рассматриваемых 
параметров субстрата ТБО, которое можно получить при сортировке. Очевидно, что 
степень извлечения отдельных сортируемых фракций при работе ЛСМ будет различ- 
ной, но для удобства расчета параметров субстрата ТБО в табл.1 и далее предполагает- 
ся, что она одинакова. 

ВЛИЯНИЕ СОРТИРОВКИ НА ТЕПЛОТЕХНИЧЕСКИЕ ХАРАКТЕРИСТИКИ 
ТБО 

Для расчета теплотехнических характеристики и других свойств мусора предполагает- 
ся, что рабочую массу ТБО составляют следующие компоненты [5]: 

С р +Н р + (У +$ р + У + А Р + \Ѵ= 100 (%) 

где С , /7°, (У, 8 , N - содержание углерода, водорода, кислорода, серы, азота (% масс.) 
в Рабочей массе ТБО; А р - зольность рабочей массы ТБО (% масс.); IV- влажность ТБО 
(% масс.). Влажность и зольность являются балластом рабочей массы ТБО. 

Состав сухой массы ТБО может быть рассчитан по формуле: 

С с + /7 е + О с +5° + іУ + А с =100 (%) ( 2) 

где С , 77 е , О , 5 е , У- содержание углерода, водорода, кислорода, серы, азота (% масс.) 
в сухой массе ТБО; А с - зольность сухой массы ТБО (% масс.). 

Следующие составляющие являются компонентами горючей массы ТБО: 

С г + Н г + о' +8 Г + Ы г - 100 (%) ( 3 ! 
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где С ,Н г ,О г , 8 Г , - содержание углерода, водорода, кислорода, серы, азота (% масс.) 

в горючей массе ТБО. 

Зная морфологический состав ТБО и содержание отдельных элементов в них, 
можно определить содержание компонентов всей массы рассматриваемых отходов. В 
дальнейшем, отношение того или иного элемента к общей массе ТБО будет помечаться 
индексом «общ», отношение элемента к отдельной фракции ТБО - индексом «/». 

Удельная низшая теплота сгорания 1 БО на рабочую массу определяется из раз- 
ницы теплоты, выделяющейся при сгорании горючих компонентов ТБО и теплоты не- 
обходимой для испарения воды, образующейся при сгорании топлива и влаги, содер- 
жащейся в ТБО. Удельная низшая теплота сгорания может быть рассчитана из форму- 
лы Менделеева [1\: 


4. 1 8 [8 1 С общ + зоо Н' общ . 26 ( О р общ 8 р бщ ) - - 6 (9 Н р бщ +ІѴ о6щ )1 кДж/кг (4), 

Г*Р ы Р глр пр 

где о6и і ’ о6и * ’ и общ > ^оби, - Общее содержание углерода, водорода, кислорода, и се- 
ры в ТБО на рабочую массу (в % масс.). ]Ѵ о6щ - общая влажность ТБО. 

Общее содержание различных химических элементов на рабочую массу ТБО 
определяется по формулам: 


= 0,01 Хсг.х, . 

І 

н об Щ = о,о\^нг 

І 

= 0 , 01 ^ 0,4 . 

і 

0,01^5,4 

/ 




где х, - содержание отдельных фракций в ТБО (в % масс); С р ,Н Р ,О р ,8 Р ,Ы Р . со- 
держание химических элементов в каждой из рассматриваемых фракций на рабочую 
массу ТБО. 

Содержание элементов в различных фракциях ТБО определяется по формулам: 


с " = с; ( 


гЛОО-ІѴ, - А р 


100 


); Щ = 


г роо-щ-л р 


I I 


100 


); 


о р = оц 


Г/ т-щ-л р 


юо 


Ц ср-сгЛОО-Щ-А,' 


( 6 ) 


М р = М‘ ( 


Г/ 100-Щ-А р 


100 


), 
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Г'Г ттГ о Г д тГ 

БД® '-д ^ і ’^і - содержание химических элементов на горючую массу от- 

дельных фракций ТБО, \Ѵі - влажность отдельных фракций ТБО; А? - зольность от- 
дельных фракций ТБО на рабочую массу. В табл. 2 приведены значения 

С, > И і > О і ,3 І , N / для различных фракций ТБО по данным [5]. 

Общая влажность и зольность на рабочую массу ТБО определяется по форму- 


лам: 


»;«,,= о, оіх>л 

/ 

44, = 0,01X44 

/ 

Общая зольность на сухую массу ТБО определяется по формуле: 

А с = А р ( і 

^общ ^обшУ+ъг. т ) 


Ш-\Ѵ 


общ 


(7) 

( 8 ) 

(9) 


Кроме 0, н , рассчитывалась также удельная низшая теплота сгорания на сухую 

массу ТБО & н ) с учетом содержания элементов в сухом веществе и отсутствия влаж- 
ности в ТБО [7]: 

- 4. 1 8 [ 8 1 С о6щ + 246 Н общ . -26 ( О общ — 8 о6щ )], кДж/кг 


Таблица 2. Содержание элементов су ,я/ ,0/ ,5/ ,У Г (% 

масс.) в горючей 

ТБО [5] 


(Ю) 


массе 


Наименование 

Содержание элементов, % масс 

фракции 

С 

Н 

О г 

У 

У" 

Бумага 

46,2 

6,2 

47,1 

0,2 

0,3 

Пищевые отходы 

53,6 

ѵг 

34,1 

0,6 

4,0 

Садовые отходы 

51,0 

6,1 

42,6 

0,1 

0,2 

Дерево 

51,0 

6,1 

42,6 

0,1 

0,2 

Резина, кожа 

77,9 

6,0 

15,1 

0,7 

0,3 

Текстиль 

56,1 

6,8 

32,2 

0,1 

4,8 

Пластмасса 

67,7 

9,3 

21,5 

0,4 

1,1 

Отсев фракций, 
менее 1 6 мм 

46,4 

6,3 

47,0 

0,3 

0 

1 
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В табл. 3 приведены значения содержания элементов в ТБО -С^ 

ПР ел д гр 

общ общ ’ 1 ' 1 общ > удельные низшие теплоты сгорания на рабочую и сухую массу ТБО; 
влажность, зольность в ТБО г.Ижевска в зависимости от степени извлечения сортируе- 
мых фракций из ТБО для различных времен года. На рис.1 приведены зависимости 

низшей теплоты сгорания на рабочую {<2 Р Н ) н сухую (0„ с ) массу ТБО от степени из- 
влечения сортируемых фракций из мусора. На рис.2 приведены зависимости влажности 
и зольности от степени извлечения сортируемых фракций из ТБО. 


Таблица 3. Свойства ТБО г. Ижевска для летнего, осеннего и зимнего периодов, в за- 
висимости от степени извлечения сортируемых фракций. Расчет сделан по 
среднему составу ТБО (см. табл.1) 


Свойства ТБО 


Степень извлечения сортируемых фракций, 

% от исходного содержания 



0 

20 

40 

60 

80 

100 

Содержание 

углерода, 

с„,% 

масс 

лето 

20,4 

20,1 

19,2 

18,6 

17,4 

15,8 

осень 

16,8 

16,3 

15,5 

14,4 

13,2 

11,6 

зима 

17,7 

17,1 

16,2 

15,2 

13,9 

12,3 

Содержание 

водорода 

н р обиг % 

масс 

лето 

2,7 

2,7 

2,6 

2,5 

2,4 

2,2 

осень 

2,1 

2,1 

2,0 

1,9 

1,8 

1,6 

зима 

2,5 

2,2 

2,1 

2,0 

1,9 

1,7 

Содержание 

кислорода 

О р - % 

лето 

14,4 

14,2 

13,6 

13,3 

12,5 

11,4 

осень 

11,4 

11,0 

10,6 

10,1 

9,4 

8,6 

масс 

зима 

12,4 

12,0 

11,4 

10,7 

9,9 

8,8 

Содержание 

серы 

•%,,,% 

масс 

лето 

0,1 

0,1 

0,2 

0,2 

0,2 

0,2 

осень 

0,1 

0,1 

0,1 

0,1 

0,1 

0,1 

зима 

0,1 

0,1 

0,1 

0,1 

0,1 

0,1 

Содержание 

азота 

лето 

0,8 

0,9 

0,9 

0,9 

0,9 

0,9 

% 

осень 

0,9 

0,9 

0,8 

0,8 

0,7 

0,6 

обиі •> /0 

масс 

зима 

0,7 

0,7 

0,7 

0,7 

0,7 

0,7 
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Влажность, 

1Т общ, %■ 

лето 

40,1 

42,6 

45,7 

48,4 

52,6 

59,2 

осень 

53,6 

56,5 

58,8 

61,7 

65,1 

69,5 

зима 

50,4 

52,8 

55,5 

58,9 

62,9 

68,1 

Зольность, 

л: бщ ,% 

масс. 

лето 

21,5 

19,4 

17,8 

16,1 

14,0 

10,3 

осень 

15,1 

13,1 

12,2 

11,0 

9,7 

8,0 

зима 

16,2 

15,1 

14,0 

12,4 

10,6 

8,3 

Зольность, 

к 6щ , % 

масс. 

лето 

35,3 

34,2 

34,0 

31,3 

29,7 

26,9 

осень 

31,9 

30,2 

29,8 

28,8 

28,1 

26,5 

зима 

32,8 

32,1 

31,6 

30,2 

28,7 

26,0 

Удельная 

теплота его- 

рания, 

3 

кДж/кг раб. 

массы 

лето 

7273,6 

7024,8 

6671,3 

6345,7 

5799,8 

5059,1 

осень 

5404,7 

5291,4 

4883,5 

4388,0 

3796,0 

3033,8 

зима 

5832,0 

5519,3 

5116,3 

4643,1 

4059,6 

3341,1 

Удельная 

теплота его- 

рания, ^ с н 

3 

кДж/кг. сух. 

массы 

лето 

14085,4 

14132,0 

14494,2 

14680,6 

15051,8 

16343,3 

осень 

15079,1 

15421,8 

15484,6 

15593,9 

15602,1 

15703,6 

зима 

14346,5 

14535,0 

14625,4 

14933,9 

15242,4 

15821,4 


Из сделанных расчетов свойств ТБО в зависимости от степени извлечения сор- 
тируемых фракций, которые могут использоваться в качестве вторичного сырья видно, 
что содержание органического вещества и общая влажность ТБО в зависимости от сте- 
пени извлечения фракций возрастают. При сортировке удельная теплота сгорания ТБО 
на рабочую массу снижается, за счет увеличения общей влажности, при этом удельная 
теплота сгорания на сухую массу ТБО несколько возрастает. 

Из сделанного анализа свойств ТБО в зависимости от степени извлечения от- 
дельных фракций, которые могут рассматриваться в качестве вторичного сырья, можно 
сделать вывод, что при сортировке ТБО ценность использования мусора в качестве то- 
плива снижается. Удельная низшая теплота сгорания ТБО на рабочую массу после 
сортировки в зависимости от времени года меньше низшей теплоты сгорания каменно- 
го угля в 6-10 раз. Одновременно увеличивается ценность использования их в качестве 
материала для получения компостов за счет увеличения содержания органического ве- 
щества и общей влажности отходов, ускоряющей процессы ферментации ТБО. За счет 
увеличения доли влажных фракций при сортировке, также увеличивается общая плот- 
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12н , кДж/кг 



■Лето 

■Осень 

■Зима 


Степень извлечения сорт, фракций, % от исх. сод. 
^ ^ II , кДж/кг 


17000 
16500 
16000 - 
15500 - 
15000 
14500 
14000 - 
13500 - 
13000 
12500 



' і 1 1 г 

0 20 40 60 80 



100 


Степень извлечения сорт, фракций, % от исх. сод. 


Рис.1. Зависимость удельной низшей теплоты сгорания на рабочую {0,^ ) и сухую 
( ^ н ) массу ТБО от степени извлечения сортируемых фракций 
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ѴУобщ, %■ 



Степень извлечения сорт, фракций, % от исх. сод. 


А р 

Л общ > 


% масс. 



Степень извлечения сорт, фракций, % от исх. сод. 

Рис.2. Приведены зависимости влажности и зольности от степени извлечения сор- 
тируемых фракций из ТБО 
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ность ТБО, что может снизить удельную площадь компостных площадок или полиго- 
нов ТБО для размещения 1 м 3 отходов. 

При уменьшении значений теплоты сгорания ТБО снижается риск самовозгора- 
ния ТБО на полигонах захоронения. Возгорание полигонов, свалок ТБО является при- 
чиной загрязнения атмосферы и других объектов окружающей среды ПХДЦ/Ф , ПХБ и 

другими токсичными веществами [8,9]. Из табл.З и рис.1 видно, что значение в 
летний период, когда наибольшая вероятность возгораний полигонов ТБО, после сор- 
тировки снижается на 20-30% и приближается к значениям 0, Р Н для несортированных 
ТБО в осенне-зимний период. 
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НАТЯЖЕНИЯ 

М. А. КОРЕПАНОВ 

Институт прикладной механики УрО РАН, Ижевск, Россия 


АННОТАЦИЯ. Предлагается метод расчета коэффициента поверхностного натя- 
жения через теплофизические свойства вещества. Предлагаемый метод позволяет 
рассчитывать коэффициент поверхностного натяжения как неполярных, так и по- 
лярных жидкостей через учет дипольного момента их молекул. Показано хорошее 
совпадение результатов расчетов коэффициента поверхностного натяжения жидко- 
стей с экспериментальными данными. 


Поверхностное натяжение входит во многие уравнения физики, физической и коллоид- 
ной химии, в частности оно определяет капиллярное давление жидкости под искрив- 
ленной поверхностью (уравнение Лапласа), давление насыщенного пара над искрив- 
ленной поверхностью (уравнение Кельвина), работу образования критического заро- 
дыша при гомогенном образовании новой фазы (уравнение Гиббса) и т.д. В то же время 
основным способом определения коэффициента поверхностного натяжения остается 
эксперимент. При этом в [1] подчеркивается, что поверхностное натяжение большинст- 
ва очень чувствительно к наличию примесей и поэтому измерения разных авторов даже 
одним методом дают неодинаковые значения сгдля одних и тех же веществ. Кроме то- 
го, для большого количества веществ коэффициент поверхностного натяжения неиз- 
вестен или известен лишь в каком-то ограниченном интервале температур, что может 
быть связано с трудностью проведения экспериментальных исследований, например, 
для криогенных жидкостей или тугоплавких соединений. 

В связи с этим актуально использование корреляционных зависимостей, позво- 
ляющих найти коэффициент поверхностного натяжения через другие свойства вещест- 
ва, известные или легко определяемые. 

Методы расчета коэффициента поверхностного натяжения рассмотрены в [2]. Так, 
Ван-дер-Ваальс, используя принцип соответственных состояний, предложил уравне- 
ние: 


(У - кТ ( ' П Р с 213 


г т У 

1 - — 


V 


( 1 ) 


У у 
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где к, п универсальные константы, Т с , Р с • - температура и давление в критической 

точке, а группа <у/Т ( ш Р ( 2ІІ может рассматриваться как безразмерное поверхностное 
натяжение. 

Еще одну корреляцию предложили Маклеод и Сагден, она связывает поверхност- 
ное натяжение с плотностями жидкости и ее пара: 

сг ' 4 = [Р\(Рж -Рп)> (2) 

где [7 3 ] - парахор, являющийся аддитивной функцией атомов и групп в молекуле и не- 
зависим от температуры. В [2] приведены таблицы со структурными составляющими 
для парахора. Корреляция Маклеода-Сагдена может быть использована для углеводо- 
родов и их производных, но неприемлема для сильно поляризованных соединений (со- 
единений с большим дипольным моментом), например, Н 2 0, НСИ и т.п 

Однако, наибольший интерес вызывают уравнения, связывающие коэффициент 
поверхностного натяжения с теплотой испарения (конденсации) вещества [3, 4]. Так в 
[э] для определения коэффициента поверхностного натяжения на поверхности раздела 
двух фаз а и /? предлагается формула 


А Н + р а ѵ а -р V + и ;(ѵ а - ѵ р ) АѴ 


где АЯ - теплота фазового перехода а ->/? , Дж/кг; ѵ - удельный объем фазы, м 3 /кг; 

и ѵ - плотность потенциальной энергии взаимодействия молекул. 

При допущении, что поверхностный слой является мономолекулярным и его 

плотность не изменена, а также р р ѵ р я КТ » р а ѵ а ,и р я 0, ѵ а /ѵ р «1 (фаза а - 

жидкость, р- газ) формула (3) в случае плоской поверхности может быть преобразова- 
на к следующему виду 

гг - АН конд~ КТ Ад 

-^ Г ’ ( 4 ) 

где т. а - координационное число, показывающее количество соседних частиц или число 
межмолекулярных связей на молекулу в конденсированной фазе; 

Ад - число недостающих соседей в поверхностном слое; г - толщина монослоя. 
Следует отметить, что формула (4) при принятых допущениях справедлива лишь 
вдали от критической точки вещества, и при ее использовании поверхностная энергия 
может иметь отрицательные значения при температурах близких к критическим. 

В то же время получить формулу для вычисления свободной поверхностной энер- 
гии аналогичную (3, 4) можно из простых рассуждений. В соответствии с «энергетиче- 
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ским» определением поверхностного натяжения свободная энергия поверхности жид- 
кости равна <7 ■ 8 , в то же время эта величина пропорциональна произведению количе- 
ства молекул в поверхностном слое жидкости на потенциальную энергию взаимодейст- 
вия молекул - Ѵ(г)-5ІЗ ,(Т) , где А, (Г) - площадь поверхности жидкости, приходя- 
щаяся на одну молекулу. 

Тогда, приравнивая эти два выражения для поверхностной энергии 

^.5 = ^), (5 ) 

можно получить выражение для коэффициента поверхностного натяжения в виде 


<6 > 

где II ( г ) - потенциальная энергия взаимодействия молекул. 

Величина потенциальной энергии взаимодействия молекул является функцией 
температуры и рода вещества и может быть описана потенциалом Леннард- Джонса [2]. 
Однако, более простым способом оценки потенциальной энергии взаимодействия мо- 
лекул, может быть, в соответствии с теорией Френкеля [5], формула 



ц А II 


хонд 


г а N. 


(7) 


где А ІІ конд - скрытая теплота испарения (конденсации), Дж/кг, которая может быть 

определена через теплоту испарения. Поскольку теплота испарения представляет собой 
сумму двух слагаемых: 

А^^ = (У / '-^)+(р Д ѵ / ’-р“ѵ«), (8) 

первое из которых является разностью удельных внутренних энергий двух фаз или 
скрытой теплотой испарения (конденсации) А , а второе - работа увеличения 

удельного объема вещества, то, скрытая теплота испарения (конденсации) NN К0Н() мо- 
жет быть определена как 

= Л Я^фѴ-рѴ) (9) 

Энергия взаимодействия молекул в жидкости II (г) в этом случае будет умень- 
шаться до нуля вслед за теплотой испарения при приближении к критической точке. 

Площадь поверхности жидкости, приходящаяся на одну молекулу А, (Т) может 
быть определена через эквимолярный радиус молекулы г, , как 
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где г, определяется как 


5,(Г) = ^ 2 


(10) 


п = 


( з V 


ѵ 4 я- р(г)л^ 


(И) 


а ѵ / ~ №І{р(Т)М - объем жидкости, приходящийся на одну молекулу, м 3 . 

Таким образом, выражение для коэффициента поверхностного натяжения может 
быть записано в виде: 


М Т \ = И Л^(г)-(рѴ-рѴ’) 

2 “ N. 


2 ' 


( 12 ) 


71 


И 


4 к р{Т )тѴ 


У 


В [4] для вычисления коэффициента поверхностного натяжения предлагается 
формула 


(7 - 


А Н 


конд 


2 71 


( 2 \і 

'ДР ' 3 


(13) 


аналогичная формулам (4, 12), отмечается качественное согласие результатов, вычис- 
ленных по формуле (13) с экспериментальными данными. 

В табл.1 приведены результаты расчетов коэффициента поверхностного натяже- 
ния по формулам (4, 12, 13) в сравнении с известными экспериментальными данными 
для различных жидкостей [1, 6]. При использовании формулы (4) толщина монослоя 
была принята равной удвоенному эквимолекулярному радиусу т = 2г г 

Из табл.1 видно, что для таких веществ как Аг, N 2 , СН 4 , С 6 Н 6 , Ыа, Ре рассчитанные 
значения коэффициента поверхностного натяжения хорошо совпадают с известными 
справочными (экспериментальными) значениями. Наибольшую погрешность дает фор- 
мула (4), что связано, во-первых, с допущениями, использованными в [3] при преобра- 
зовании формулы (3) в (4), а во-вторых, с тем, что поверхностный монослой может 
иметь структуру, отличающуюся от структуры всей жидкости. В результаты расчетов 
по формуле (13) вносится погрешность, связанная с тем, что в ней не учитывается ко- 
ординационное число жидкости (число соседних атомов), которое для различных ве- 
ществ может составлять от 4 до 8 и более. 

Однако, в [2, 4] отмечается, что существующие методы расчета поверхностной 
энергии (коэффициента поверхностного натяжения) дают большую погрешность для 
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Таблица 1. Сравнение методов расчета коэффициента поверхностного натяжения 


Вещество 

Т, К 

р, кг/м - ' 

АН, 


а, Дж/м л 





кДж/кг 

форму- 
ла 4 

форму- 
ла 12 

форму- 
ла 13 

справ. 

Аргон, Аг 

87 

1407.5 

162.7 

0.0152 

0.0101 

0.0132 

0.0113 

Азот, N 2 

77 

808.0 

199.5 

0.0114 

0.0076 

0.0099 

0.00885 


90 

746.3 

173.8 

0.0085 

0.0060 

0.0082 

0.00616 

Аммиак, 

Ш 3 

240 

681.7 

1369.8 

0.0917 

0.0613 

0.0515 

0.031 

Вода, Н 2 0 

373 

958.13 

2257.2 

0.1774 

0.1148 

0.1084 

0.0589 

Метан, СН4 

109 

466.0 

509.88 

0.0168 

0.0112 

0.0146 

0.0137 

Бензол, 

с 6 н 6 

355.2 

825.7 

393.94 

0.0328 

0.0219 

0.0279 

0.0212 

283 

877.2 

512.11 

0.0462 

0.0308 

0.0378 

0.03024 

Натрий, Ыа 

1159 

739.0 

3919.1 

0.1992 

0.1328 

0.1717 

0.1278 

Железо, Ре 

1811 

7024 

7502.0 

2.4828 

1.659 

1.982 

1.850 

Хлорметан, 

СНзСІ 

250 

915.9 

428.42 

0.0330 

0.0221 

0.0281 

0.0120 


полярных жидкостей, т.е. жидкостей, молекулы которых обладают большим диполь- 
ным моментом. Для подобных жидкостей коэффициенты поверхностного натяжения, 
вычисленные по формулам (4, 12, 13) имеют существенно завышенные значения по 
сравнению с экспериментальными данными, что видно из табл.1 (Н 2 0, РІН 3 , СН 3 С1). 

Объяснение этому можно найти в природе вандерваальсовых сил, действующих 
между молекулами на далеких расстояниях: 

1) . Электрические силы между молекулами, обладающими постоянными электри- 
ческими диполями. Взаимодействие этого типа известно под именем ориентационного 
эффекта. 

2) . Силы, вызываемые взаимодействием индуцированных диполей - индукцион- 
ный эффект. 

3) . Дисперсионный эффект, присущий всем газам. 

Все эти три силы являются силами притяжения и в первом приближении могут 
рассматриваться как аддитивно складывающиеся центральные силы [5]. В табл. 2 при- 
ведены все три составляющие молекулярных сил для некоторых веществ. Из нее видно, 
что чем больше дипольный момент молекулы, тем больший вклад в энергию молеку- 
лярного взаимодействия дает ориентационный эффект. А у молекул с очень малым ди- 
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Таблица 2. Три составляющие Вандерваальсовой энергии взаимодействия молекул 


Вещество 

Дипольный 
момент, Д 

Ѵ ор 'ІО 60 , эрг 

Ѵинд 'ІО 60 , эрг 

Ѵдисп 'ІО 60 , Эрг 

СО 

ОЛ 

0,0034 

0,057 

67,5 

Ш 

0,5 

0,35 

1,68 

3,82 

НС1 

1,1 

18,6 

5,4 

105 

N43 

1,5 

84 

10 

93 

н 2 о 

1,8 

190 

10 

47 


польным моментом, например СО, или неполярных молекул или атомарных веществ, 
энергия взаимодействия определяется исключительно дисперсионным эффектом. 

С другой стороны, в соответствии с основными принципами термодинамики, сис- 
тема, в данном случае - жидкость, стремится к состоянию с минимальной энергией, т.е. 
жидкость должна минимизировать свободную поверхностную энергию Р=о5 [9]. Этим 
можно объяснить и то, что капля жидкости принимает шарообразную форму - вещест- 
во стремится минимизировать площади поверхности, соответствующую своему объе- 
му. Кроме того, можно предположить, что в полярной жидкости молекулы 
в поверхностном слое ориентируются таким образом, чтобы уменьшить свободную по- 
верхностную энергию за счет дипольного взаимодействия, т.е. ориентационного эф- 
фекта. 

Для учета ориентационного эффекта в полярных жидкостях в [10, 11] использует- 
ся потенциал Штокмайера, который может быть выражен через потенциал Леннард- 
Джонса и энергию дипольного взаимодействия молекул [10]. Использование потенциа- 
ла Леннард-Джонса затруднено необходимостью знать глубину потенциальной ямы и 
«эффективное» межмолекулярное расстояние вещества. В то же время представляется 
возможным вместо потенциала Леннард-Джонса использовать значение энергии взаи- 
модействия молекул определяемое выражениями (7, 9). 

Таким образом, для расчета коэффициента поверхностного натяжения полярных 
жидкостей энергию взаимодействия молекул, вычисленную, например, по формуле (7) 

необходимо уменьшить на энергию дипольного взаимодействия молекул, значение ко- 
торой равно: 

и = р 2 _ 3 ЫЫ 

дш 4/гб- 0 г 3 4 тге 0 г 5 ’ ^ 14) 

где р — дипольный момент молекулы; 
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г - расстояние между молекулами, которое может быть принято равным удвоенно- 
му эквимолярному радиусу 2г, (11). 

Максимальное значение энергии дипольного взаимодействия достигается в случае 
встречного расположения диполей “+” к 


у - Р 

и дип 3 

7ІБ { / 


( 15 ) 


Таким образом, для расчета коэффициента поверхностного натяжения жидкостей 
возможно использование следующей формулы 


г(Г) = 


И А Я„„Дг) + р“ѵ“-рѴ 


N 


-к- 


7ГЕг. 


(2г,У 


(16) 


где к - коэффициент, учитывающий тепловое движение молекул и удовлетворяющий 
условию Ііш к = 0 . 

Т-^Т 

* г л кр 

В табл.З. приведены значения коэффициента поверхностного натяжения для неко- 
торых полярных жидкостей в широком температурном диапазоне. 

При расчете коэффициента поверхностного натяжения по формуле (16) коэффи- 
циент к был принят постоянным во всем температурном диапазоне и равным 1/3 для 
большинства веществ, кроме гидразина и оксида свинца. 

Из формулы (15) следует, что величина поправки на дипольное взаимодействие зависит 
от дипольного момента молекул и межмолекулярного расстояния. Из таблицы видно, 
что для таких веществ как хлорбензол и анилин поправка дает уточнение коэффициен- 
та поверхностного натяжения примерно на 1 0% по сравнению с результатами по фор- 
муле (12). В то же время для аммиака и воды - на 50% и более при сравнимом диполь- 
ном моменте молекулы, что объясняется в первую очередь малостью молекул аммиака 
и воды по сравнению с хлорбензолом или анилином. Также можно отметить, что для 
глицерина и нитробензола поправка дает уточнение на 20-30%, это при сопоставимых 
размерах молекул последних с хлорбензолом и анилином объясняется вдвое большим 
дипольным моментом. 

При этом следует отметить, что в работе [10], делается вывод о практически пре- 
небрежимом влиянии ориентационного эффекта на поверхностное натяжение, что про- 
тиворечит полученным результатам. 

Сравнение результатов расчетов с экспериментальными данными [1,6] показыва- 
ет хорошее совпадение для различных веществ в широком диапазоне температур. 
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Таблица 3. Коэффициент поверхностного натяжения 


Вещество 

Т, К 

Р, 

кг/м 3 

АН, 

кДж/к 

г 

Диполь- 
ный мо- 
мент, Д 

к 

а, Дж/м 2 

форм. 

3.20 

форм. 

3.23 

справ. 

Аммиак, ЫНз 

240 

681.7 

1369.8 

1.50 

1/3 

0.06130 

0.03510 

0.03100 

284 

625.7 

1170.1 

1.50 

1/3 

0.04900 

0.02640 

0.02340 

Гидразин, 

ы 2 н 4 

293 

1004.0 

1684.2 

1.80 

1/4 

0.08361 

0.06484 

0.06670 

386.5 

970.0 

1397.0 

1.80 

1/4 

0.06590 

0.04817 

0.04900 

Вода, Н 2 0 

273 

999.87 

2501.0 

1.80 

1/3 

0.13859 

0.07368 

0.07565 

373 

958.13 

2257.2 

1.80 

1/3 

0.11480 

0.05790 

0.05890 

473 

864.9 

1939.0 

1.80 

1/3 

0.09234 

0.04137 

0.03768 

573 

712.2 

1403.0 

1.80 

1/3 

0.05721 

0.02034 

0.01437 

623 

574.5 

722.1 

1.80 

1/3 

0.02912 

0.00335 

0.00368 

Хлорметан, 

СН 3 С1 

250 

915.9 

428.4 

1.90 

1/3 

0.02210 

0.01089 

0.01200 

Метанол, 

СНзОН 

273 

810.0 

1210.0 

1.70 

1/3 

0.05133 

0.03571 

0.02450 

373 

714.0 

1030.0 

1.70 

1/3 

0.03867 

0.02602 

0.01570 

Глицерин, 

С 3 Н 8 0з 

293 

1260.0 

968.8 

3.00 

1/3 

0.08226 

0.06479 

0.06300 

423 

1155.0 

1313.6 

3.00 

1/3 

0.05671 

0.04160 

0.04880 

Этиловый 

эфир, 

С 2 Н 5 ОС 2 Н 5 

303 

701.9 

353.8 

1.30 

1/3 

0.01729 

0.01551 

0.01580 

423 

517.9 

215.6 

1.30 

1/3 

0.00828 

0.00721 

0.00310 

Хлорбензол, 

С 6 Н 6 С1 

405 

980.0 

324.9 

1.60 

1/3 

0.02295 

0.02061 

0.02070 

Нитробензол, 

СбН5ГЮ 2 

293 

1203.0 

678.0 

3.10 

1/3 

0.06047 

0.04982 

0.04390 

483 

1100.0 

374.0 

3.10 

1/3 

0.02956 

0.02038 

0.02360 

Анилин, 

с 6 н 5 ш 2 

457 

875.0 

448.0 

1.60 

1/3 

0.02757 

0.02491 

0.02400 

Оксид угле- 
рода, СО 

81.5 

789 

215.8 

0.10 

1/3 

0.00811 

0.00797 

0.00800 

Оксид желе- 
за, РеО 

1647 

4600 

6742.2 

4.70 

1/3 

1.22979 

0.66822 

0.58500 

Оксид свин- 
ца, РЬО 

1159 

7600 

1055.4 

4.67 

1/2 

0.40399 

_ 

0.11367 

0.13200 
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РАСЧЕТ ТЕРМОДИНАМИЧЕСКИХ СВОЙСТВ МЕТАСТАБИЛЬНЫХ 
КОНДЕНСИРОВАННЫХ КЛАСТЕРОВ ІУВ2ГѴВ5 И СОСТАВА ЖИДКОГО 
НИОБИЯ С УЧЕТОМ ИХ СУЩЕСТВОВАНИЯ 

Г.К. МОИСЕЕВ 

Российская академия наук, Уральское отделение, ГУ Институт металлургии, 
Россия, Екатеринбург 


АННОТАЦИЯ. С использованием данных молекулярной динамики и известных 
расчетных методов оценены термодинамические свойства метастабильных конден- 
сированных кластеров - N65; с учетом атомов и кластеров методами термоди- 
намического моделирования рассчитаны составы жидкого ниобия при 2800-5300 К; 
Т и АН кипения. 


Изучение состава расплавов индивидуальных ЩМ, А1, 8і и их смесей показало, что в 
них, кроме атомов, существуют "малые" метастабильные кластеры, в том числе самоас- 
социаты [1-3]. Известны термодинамические данные для моноатомных конденсирован- 
ной и газовой фаз ниобия [4]; равновесных характеристик метастабильных самоассо- 
циатов ниобия нами не обнаружено. 

Целями настоящей работы являлось: - оценка термохимических свойств и тер- 
модинамических функций некоторых метастабильных "малых" кластеров ниобия; - 
расчет состава жидкого ниобия в широком интервале температур с учетом атомов и 
малых кластеров; Т и АН кипения модельного расплава-раствора ниобия. 

ТЕРМОДИНАМИЧЕСКИЕ СВОЙСТВА И ФУНКЦИИ 

Для определения величин АН“ ;1( (СЭО) конденсированных кластеров использовали 
средние энергии связей между атомами в изолированных кластерах ниобия (|Е Ь |, 
эВ/атом), полученные для к*ЫЬ п с п = 2+23 методами молекулярной динамики [5]. Рас- 
четы проводили по соотношению [6] 

АН» 98 (к*ЫЬ п ) « 59401 -п-|Еь| (№>„), Дж/моль, (1) 

и получили зависимость 
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ДН° 98 (к*М) п ) = -33.984 + 30. 179-п, кДж/моль, (2) 

со средним отклонением ±1.55%. Анализ (2) показал, что СЭО кластеров могут изме- 
няться в пределах 26.4 (п = 2) -705 (п=23) кДж/моль. Для последующих исследований 
использовали СЭО (к*ЫЬ п ) с п = 2-5, т.е. минимальные по величине энергетики обра- 
зования самоассоциаты. 

Для расчетов СЭО конденсированных кластеров использовали также соотноше- 
ние [6] 


ДН« 98 (к*ИЪ„) * 59401 -п-ІЕьКСп)- ^°^ Ь) , Дж/моль, 


( 3 ) 


где |Еь|(С п ) взято по данным [7], а величины электроотрицательностей (ЭО) по [8]. 

Исходные величины |Еь|( БИу), СЭО кластеров с п = 2-5, рассчитанные по (1) и 
(3) приведены в табл. 1. Среднеарифметические величины СЭО приняты за достовер- 
ные. Остальные термохимические характеристики рассчитаны с использованием из- 
вестных данных для к*ЫЬі [4] и апробированных методов расчета [9, 10]. 

Принятые за достоверные характеристики и температурные зависимости приве- 
денной энергии Гиббса даны в табл. 2. 

Сведений о существовании летучих самоассоциатов ниобия нам не обнаружено. 
Поэтому компьютерные эксперименты выполнены с привлечением только конденсиро- 
ванных "малых" кластеров ниобия. 


ТЕРМОДИНАМИЧЕСКОЕ МОДЕЛИРОВАНИЕ (ТМ) 

Выполнено с использованием методологии ТМ [11], пакета АСТРА. 4 с БД 
АСТРА.ВА8 [12]. Исходное рабочее тело "99 масс % N6 +1 масс % Аг"; расчеты вы- 
полняли при 2800 - 5300 К с шагом 100° (в области 5100-5300 К с шагом 10°), общем 
давлении 9.8066Т0' 2 МПа. Для описания расплава использовали модель идеального 
раствора продуктов взаимодействия (ИРПВ) [11], апробированную при исследовании 
расплавов щелочных металлов [1]. В составе ИРПВ учитывали присутствие 
[ІГЪі]-[№> 5]; в составе газовой фазы {N6]}, {Аг}, {Аг +1 } и е-газ. 


Таблица 1. Величины |Еь| для к*(74Ь2 - N65) и расчетные значения их СЭО 


п 

|Еь|, 

эВ/атом, [5] 

СЭО, кДж/моль 

(1) 

(3) 

среднее 

2 

-2.22 

26.374 

24.71 

25.542 ±0.832 

3 

-3.03 

53.995 

58.16 

56.077 ±2.082 

4 

-3.74 

88.864 

80.6 

84.732 ±4.132 

5 

-3.96 

117.614 

106.45 

1 12.032±5.582 
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Таблица2. А. Принятые за достоверные термохимические характеристики метастабильных 
конденсированных самоассоциатов N62 - N65 


РАСЧЕТ ТЕРМОДИНАМИЧЕСКИХ СВОЙСТВ МЕТАСТАБИЛЬНЫХ КОНДЕНСИРОВАННЫХ КЛА- 
СТЕРОВ ИВз-ИВз И СОСТАВА ЖИДКОГО НИОБИЯ С УЧЕТОМ ИХ СУЩЕСТВОВАНИЯ 
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РЕЗУЛЬТАТЫ И ОБСУЖДЕНИЕ 

Состав модельного расплава показан на рисунке. Жидкий ниобий содержит атомы и 
"малые" кластеры. При всех температурах наблюдается иерархия их содержаний: 
х [№>]]> х[ЛЬ 2 ]> х[№>з]> х[№> 4 ]> х[ЫЬ 5 ], т.е. вероятность образования и существования 
названных частиц уменьшается в той же последовательности. Особый интерес вызыва- 
ет изменение содержания частиц с ростом температуры, где выявляются заслуживаю- 
щие внимания особенности. Рассмотрим "поведение" кластеров. Рост температуры со- 



Рис. Состав модельного расплава ниобия согласно ТМ. 

1 - х [ЫЬі], 2 - х РМЬ 2 ], 3 - х [№> 3 ], 4 - х [N 64 ], 5 - х [ЛЬ 5 ] 
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РАСЧЕТ ТЕРМОДИНАМИЧЕСКИХ СВОЙСТВ МЕТАСТАБИЛЬНЫХ КОНДЕНСИРОВАННЫХ КЛА- 
СТЕРОВ ^ 2 -ЫВ 5 И СОСТАВА ЖИДКОГО НИОБИЯ С УЧЕТОМ ИХ СУЩЕСТВОВАНИЯ 


провождается увеличением х[ЫЬ 5 ] и х[Е1Ъ 4 ] во всем интервале температур до 5200 К; 
х[№>з] до -4600 К и х[№ 2 ] до -3600 К. затем при 4600-5200 К х[ЫЬз] уменьшается на 
-0.1 мол %; при 3600-5200 К х[ЫЬ 2 ] уменьшается на -1.5 мол %. 

Атомная составляющая расплава, - х[№>і], -2800 - 3800 К (4000 К) уменьшается 
от 75.55 до - 73.0 мол. %, что объясняется ростом содержания кластеров. При 4000 - 
5200 К х[ЫЬі] растет от 73 % до 74 %. Увеличение содержания [№>і], по-видимому, 
связано с диссоциацией кластеров [№>2] и [N63] по реакциям: 

[НЬ 2 ] = 2[ЫЬ,]. (4) 

[НЬ 3 ] = ЗрЧЪ,]. (5) 

Из табл. 1 следует, что |Е Ь | [ИЬ 5 ] : |Е Ь | [ЫЬ 4 ] : |Е Ь | [№> 3 ] : |Е Ь | [ЫЬ 2 ] = 

1.784:1.685:1.365:1. Поскольку исследования выполнены при высоких температурах, то 
обменная энергия (или энергия взаимодействия) между частицами в расплаве ниобия 
является значительной и при определенных температурах будет превышать средние 
энергии связей между атомами в кластерах. Это является наиболее вероятной причиной 
реализации в первую очередь превращений (4) и (5). Сопоставление температур начала 
диссоциации кластеров (начала уменьшения х [ЗМЬ 2 ] и х [N63] по рис.). Т § [ЫЬ 3 ] / Т 8 [N62] 
= 4600 / 3600 = 1.278 коррелирует с отношением |Е Ь | [N63] / |Е Ь | [№>2] = 1.365. Поэтому 
в первую очередь наблюдается диссоциация по схеме (4) и по мере последующего рос- 
та температуры - по схеме (5) (кривые 2 и 3 рисунка). Поскольку "прочность" класте- 
ров [ЫЬ 4 ] и [N65] значительно выше, чем [N62] и [N63], то в изученном интервале 
температур пороговые" величины обменной энергии в расплаве не достигаются, и не 
наблюдается уменьшения содержания этих частиц. 

Термодинамическая возможность превращений (4) и (5) подтверждается оцен- 
кой величин их ДОа (табл. 3). 

Расчетная температура кипения равна 5190 ± 10° К, ДН кипения « 662 кДж/моль. 
По данным [8, 13] Т кипения равна 5115 и 5136 К, соответственно, и ДН кипения = 
683.7 кДж/моль [8]. 

Таблица 3. Оценка величин ДО реакций (4) и (5) 


Реакция 

ДО , Дж/моль [Мц], при Т, К 

298 

3000 

3600 

5000 

[ЫЬ 2 ] = 2 [N61] 

38,2 

-192,5 

-2000 

-9240 

[N63] = зрчгь,] 

39,0 

-198,0 

-6760 

-10560 
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ЗАКЛЮЧЕНИЕ 

По-видимому, впервые оценены термодинамические свойства и функции конденсиро- 
ванных метастабильных кластеров N 62 - 5 . Рассчитан состав модельного расплава нио- 
бия, содержащий атомы и кластеры, при 2800-5300 К, Р = 908066- 10“ МПа в исходной 
среде Аг; Т и ДН кипения расплава. Обсуждены зависимости содержания частиц в рас- 
плаве от температуры. Расчетные характеристики кипения хорошо согласуются со 
справочными данными. 

Работа выполнена при финансовой поддержке программы исследований Прези- 
диума РАН "Фундаментальные проблемы физики и химии наноразмерных систем и на- 
номатериалов"; РФФИ (проект № 04-03-33109). 
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АННОТАЦИЯ. В статье рассмотрена полученная с применением метода молеку- 
лярной динамики модель формирования острия зондирующей иглы С ГМ. Приведе 
ны результаты исследования влияния силы тяжести нижней части заготовки иглы 
на характер разрыва ее «шейки» с целью получения атомарных микровыступов в 
области разрыва. 


ВВЕДЕНИЕ 

Современная туннельная микроскопия является одним из немногих методов исследо- 
вания поверхности, позволяющим производить неразрушающий контроль ее парамет- 
ров с атомным разрешением. Для получения предельной разрешающей способности 
сканирующего туннельного микроскопа (СТМ) его зондирующая игла (ЗИ) должна 
иметь атомарную остроту. Известно, что при изготовлении ЗИ (например, методами 
химического или электрохимического перетравливания проволоки-заготовки) в месте 
разрыва «шейки» заготовки существует вероятность образования микровыступов, 
имеющих размеры -20-200 А. Поэтому актуальной является задача создания методик, 
позволяющих заострять подобные микровыступы до атомарно острого состояния. 

В настоящее время наиболее перспективными методами заострения микровы- 
ступов являются полевые методы [1,2]. Они обладают тем преимуществом, что с их 
помощью можно изготавливать атомарно острые ЗИ непосредственно в условиях СТМ 
(іп 5 ІШ). Это означает, что для контроля остроты ЗИ не применяются методики, исполь- 
зование которых сопряжено со значительными экспериментальными трудностями и 
затратами (например, обеспечение высокого вакуума, низких температур и т.д). Недос- 
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татком такого подхода является то, что полевое заострение ЗИ происходит «вслепую», 
поскольку экспериментатор не обладает точными данными о наличии и свойствах мик- 
ровыступов, формирующих атомное пространственное разрешение СТМ. 

Поэтому для создания и улучшения методик атомарного заострения ЗИ іп зіш 
необходима информация о характере образования и точной топологии микровыступов, 
образующихся на «изломе» заготовки ЗИ. Кроме того, необходимыми являются реко- 
мендации для усовершенствования методов изготовления ЗИ, имеющих микровыступы 
оптимальной для последующего заострения формы. 

При изготовлении ЗИ методами химического или электрохимического травления 
в растворе электролита реакция проходит однородно по всей поверхности погруженной 
части проволочной заготовки, за исключением мениска на поверхности электролита, 
где реакция протекает быстрее. В итоге в этой области заготовки образуется «шейка», 
диаметр которой постепенно уменьшается, и в результате нижняя часть заготовки от- 
рывается. Задачей данного исследования является изучение влияния силы тяжести на 
характер разрыва «шейки» заготовки ЗИ СТМ, а также на формирование микровысту- 
пов на месте излома. 

МОДЕЛЬ ЗИ СТМ 


В качестве метода исследования был использован метод молекулярной динамики. Его 
применение обусловлено уникальностью, позволяющей исследовать поведение слож- 
ных систем, состоящих из сотен тысяч атомов. С другой стороны, моделирование по- 
добных кластеров имеет недостатки, выражающиеся: во-первых, в применении методов 
классической физики (а не квантовой механики, что потребовало бы огромных вычис- 
лительных мощностей), а во-вторых, в использовании исключительно эмпирических 
параметров. Тем не менее, молекулярно-динамические исследования имеют сущест- 
венное значение при исследовании свойств кластерных материалов [3,4]. 

В молекулярно-динамических расчетах будущие положения и скорости атомов 
определяются их текущими скоростями и положениями. При этом сила, действующая 
на і -й атом, определяется градиентом потенциальной энергии: 



( 1 ) 


где V - функция потенциальной энергии; г ( - положение і -го атома. Поведение сис- 
темы из N атомов определяется при помощи системы дифференциальных уравнений: 



дѴ 


(2) 
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где т і - масса 1-го атома; І — Х...Ы . 

Для исследования процессов разрыва «шейки» заготовки ЗИ необходимы изуче- 
ние и анализ деформации и разрыва кластера, моделирующего область разрыва под 
действием силы тяжести. В качестве объекта моделирования использовался кластер из 
19949 атомов вольфрама, расположенных в его кристаллической решетке (рис. 1). 
Предполагается что, заготовка состоит из «шейки», а также верхней и нижней массив- 
ных частей, схематично представленных на рис. 2. При этом масса «шейки» является 
пренебрежимо малой по сравнению с массой нижней части. Поэтому вес М§ нижнего 
основания является единственным параметром (в рассматриваемых условиях), влияю- 
щим на процесс разрыва «шейки» заготовки. Это предположение является вполне 
обоснованным, так как при электрохимическом травлении (например, в «петлевом» от- 
рицательном электроде, заполненном электролитом [5]) наиболее интенсивное утонь- 
шение проволоки-заготовки происходит именно в области ее «шейки». Поэтому сила 
тяжести, обусловленная весом нижней части заготовки, при моделировании распреде- 
ляется на атомы «шейки», что обеспечивает ее деформацию и, в конечном итоге, раз- 
рыв в наиболее тонком месте. 



Рис Л. Кластер «шейки» заготовки из 19949 атомов вольфрама 
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Т 


Рис. 2. Схематичное изображение заготовки зондирующего острия СТМ в момент 
разрыва 

Известно, что малоразмерные металлические объекты (кластеры) обладают 
структурой, несколько отличающейся от структуры кристаллических металлов. При 
этом существенность отличий зависит от размеров кластера[6,7]. Атомы «шейки», из- 
начально располагавшиеся в ячейках кристаллической решетки вольфрама, смеща- 
ются из них, занимая новые положения, отвечающие принципу наименьшего дейст- 
вия. Поэтому исследования разрыва «шейки» заготовки должны предваряться распо- 
ложением ее атомов в равновесные положения. 

РЕЗУЛЬТАТЫ И ОБСУЖДЕНИЕ 

Для моделирования разрыва «шейки» заготовки ЗИ СТМ методом молекулярной дина- 
мики использовался программный пакет ИАМО, созданный в Иллинойском универси- 
тете (США) . ЫАМО хорошо зарекомендовал себя при расчетах динамики кластеров и 
ультрадисперсных частиц, состоящих из большого числа атомов [8]. Такие опциональ- 
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ные возможности как, например, расчет молекулярной динамики при наличии термо- 
стата, возможность оптимизации геометрии кластера путем минимизации энергии, оп- 
ределение различными способами действующих на атомы сил, а также бесплатное рас- 
пространение (Тгее\ѵаге) и открытость кода делают пакет ЫАМО оптимальным инстру- 
ментом для исследования разрыва «шейки» заготовки ЗИ СТМ. 

Действие силы тяжести на «шейку» заготовки моделировалось следующим об- 
разом. 

1. Производилась оптимизация геометрических положений атомов «шейки» пу- 
тем минимизации энергии. 

2. При помощи экспериментальных данных рассчитывался вес М§ нижней час- 
ти заготовки (рис. 2), которым (путем перераспределения на отдельные атомы) нагру- 
жалась «шейка» (при этом атомы ее верхней части считались неподвижно закреплен- 
ными). 

3. Выполнялось молекулярно-динамическое моделирование разрыва «шейки» 
заготовки. 

При расчетах использовались следующие эмпирические данные: длина нижней 
части заготовки I - 0.2 см; радиус цилиндрической нижней части заготовки 

К = 0.002 см; плотность вольфрама р = 19.3 г/ см 3 . При этом сила, действующая на 
атом «шейки», вычислялась согласно выражения: 


р = рлЯ^1 

N 


( 3 ) 


где іѴ- число «незакрепленных» атомов «шейки»; ^-ускорение свободного падения. 
Потенциал межатомного взаимодействия Деннарда- Джонса определялся как: 




(4) 





где г ( - расстояние между атомами типа і, на котором потенциальная энергия их взаи- 
модействия минимальна; 8 1 - глубина потенциальной ямы для двух атомов типа і. При 
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этом значения эмпирических параметров для вольфрама принимались равными: 
г = 1 .650 А; 8 1 =0.150 ккал/моль [8]. 

Сила, обусловленная весом нижней части заготовки и вычисленная согласно вы- 
ражения (3), оказалась равной 1.75-10 'ккалДмоль-А). Будучи приложенной к атомам 
«шейки», она оказалась недостаточной для ее разрыва. Более того, численные экспе- 
рименты показали, что для разрыва «шейки» необходимы силы, примерно на порядок 
превышающие вес нижней части заготовки ЗИ СТМ (возникающие, например, при уд- 
линении и утяжелении нижней части). При этом микротопография области такого раз- 
рыва заготовки ЗИ представлена на рис. 3. 

Таким образом, можно сделать вывод, что разрыв «шейки» заготовки ЗИ СТМ 
при ее изготовлении методом электрохимического травления невозможно объяснить 
лишь действием силы, обусловленной весом нижней части заготовки. Природа этого 
процесса более сложна, так как требует учета сил значительно больших, чем вес этой 
нижней части. Такие силы действительно присутствуют в момент разрыва. Они обу- 
словлены колебаниями нецентрированной нижней части заготовки, а также ее колеба- 
ниями под действием ударов молекул раствора электролита, находящихся в постоян- 



Рис. 3. Микротопография кончика ЗИ СТМ, образованной в результате разрыва 
«шейки» заготовки (вид снизу) 
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ном тепловом движении. Эти колебания становятся весьма значительными при силь- 
ном (-100 А) утонынении «шейки» заготовки. В целом, процесс формирования острия 
ЗИ можно считать следствием воздействия таких сил и силы тяжести нижней части за- 
готовки, приводящих к совместному процессу излома и разрыва «шейки» заготовки. 
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АННОТАЦИЯ. Приведены результаты исследований воздействия постоянного 
магнитного поля на флюид нефть плюс вода с целью построения эффективных маг- 
нитных устройств. Показано, что применение магнитной обработки позволяет су- 
щественно уменьшить образование отложений на стенках насосно-компрессорных 
труб, и увеличить срок работы скважин. 


Совершенствование техники и технологии добычи нефти повышенной вязкости являет- 
ся одной из наиболее актуальных задач нефтедобывающей отрасли. Сложность в реше- 
нии этой проблемы заключается в многообразии разрабатываемых объектов и место- 
рождений, характеризующихся высоким содержанием асфальтосмолопарафиновых 
(АСП) компонентов в добываемом продукте. При подъеме такого продукта и его охла- 
ждении происходит интенсивное отложение этих асфальтосмолопарафиновых состав- 
ляющих, прежде всего на стенках насосно-компрессорных труб (НКТ). Для восстанов- 
ления эффективности таких скважин необходима механическая чистка или горячая 
промывка. ѵ 

Еще в 1960 г. показали, что магнитная обработка позволяет уменьшить отложе- 
ния не только солей жесткости, но и существенно уменьшить образование асфальтас- 
молопарафиновых отложений (АСГІО) [1]. Позже было показано, что магнитная обра- 
ботка жидкостей существенно меняет их коррозионную активность (см. например[2]) 

К концу прошлого века научно-производственные предприятия ряда городов 
России начали выпускать устройства для магнитной обработки жидкости, в том числе и 
нефти при ее добыче и транспортировке [3]. К настоящему времени вопросам воздей- 
ствия магнитного поля на добываемую нефть посвящено достаточно большое количе- 
ство научно-технических публикаций [4-8]. Преимущественная часть их касается прак- 
тического использования магнитных систем. При этом отмечаются высокие показатели 
данного способа обработки нефти. Эффективность находится на уровне лучших тради- 
ционных средств при значительно меньшей стоимости (в десятки раз). Устройства спо- 
собны работать длительный период, совпадающий со сроком службы НКТ, и без затрат 
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в процессе эксплуатации. Особо следует отметить экологическую чистоту магнитного 
способа обработки нефти. 

Механизмы же воздействия магнитного поля на добываемый флюид изучены 
крайне слабо. 

В отделе магнетизма твердых тел НИИ физики и прикладной математики 
Уральского госуниверситета совместно с ЦНТ «ЛАНТАН», проводятся работы по изу- 
чению механизмов воздействия постоянного магнитного поля на флюид нефть плюс 
вода с целью построения эффективных магнитных устройств для различных типов 
нефти. 

Напомним, что основными компонентами отложений на внутренней поверхно- 
сти НКТ являются сера, парафин, асфальтены и смолы. В пласте и призабойной зоне 
скважины при достаточно высокой температуре они как бы растворены в основных уг- 
леводородных составляющих нефти и равномерно распределены по объему нефтеводя- 
ной эмульсии. При подъеме нефти по НКТ уменьшается ее температура и на каком-то 
уровне она опускается до температуры кристаллизации этих компонентов. Естественно, 
что, прежде всего, температура кристаллизации достигается у стенок НКТ, поэтому на 
внутренней поверхности НКТ и появляются зародыши кристаллизации. 

Прежде чем говорить о влиянии магнитного поля на нефть, попытаемся опреде- 
лить, с какими в магнитном отношении компонентами продукта мы имеем дело. Изме- 
рение магнитной восприимчивости х ряда нефтей Удмуртии и Когалыма показывают ее 
диамагнетизм, а величины отрицательной восприимчивости лежат в интервале (-б - ^- 

9) 10 тогда как отдельные нефти Татарстана и Перми имеют восприимчивость — 
0,95- ІО' 6 . 

В книге «Магнитные свойства нефтей» [9] приводятся значения магнитной вос- 
приимчивости нефтей от -0,75- ІО" 6 до -0,83Т0" 6 . Авторы утверждают, что ими « иссле- 
дованы магнитные свойства нефтей большого количества месторождений различных 
по возрасту нефтевмещающих отложений и приуроченности к геотектоническим 
структурам: Волго-Уральской, Тимано-Печорской, Предкавказья и Северного Сахали- 
на». 

Следует отметить, что при многообещающей заявке - «исследованы магнитные 
свойства нефтей (более 1000 образцов)» - на самом деле измерялась лишь магнитная 
восприимчивость с сомнительной точностью, т.к. разница в значениях восприимчиво- 
сти составляет менее 10%, несмотря на разный возраст и провинции, а следовательно 
различное содержание диамагнитных, парамагнитных и ферромагнитных компонент, 
определяющих величину и знак восприимчивости флюида. 

Так например, ловушечная нефть Ельниковской УПН имеет положительное зна- 
чение восприимчивости 1510 6 . После ультразвуковой обработки отделилось -12% во- 
ды с отрицательной восприимчивостью -0,53Т0" 6 , а восприимчивость оставшейся неф- 
ти увеличилась до 16,3 ТО' 6 . 
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Для более глубокого изучения магнетизма добываемой нефти, ее образцы были 
исследованы на сверхчувствительном СКВИД - магнитометре ИФМ УрО РАН. На 
рис. 1 представлены зависимости удельной намагниченности двух образцов нефти 
Гремихинского месторождения от напряженности магнитного поля. Они имеют вид, 
присущий диамагнетику, с малым содержанием ферромагнитной составляющей. Урав- 
нения описывающее эти зависимости представлены в квадратной рамке на рисунке. 

Коэффициент А уравнения соответствует отрезку, отсекаемому прямой на оси 

намагниченности и определяет ферромагнетизм нефти. Среднее значение А«2,5-10' 5 
ети/§. 

Для оценки количественного содержания примесей железа в нефти, воспользу- 
емся данными работы [8]. Из них следует, что примеси сформированы в основном фер- 
ромагнитными микрокристаллами оксидов (гидроокислов железа). И они присутствуют 
в виде агрегатов массой порядке ІО" 14 г. Для расчетов примем, что ферромагнетизм 



Рис. 1 . Зависимость удельной намагниченности нефти Гремихинского месторож- 
дения от напряженности магнитного поля 
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нефти определяется магнетитом РезСЩ который имеет удельную намагниченность 
Мрезо 4 = 92 еши/§. (Ферромагнетизм, Р.Бозорт стр. 1 93. М. 1956г.) 

Тогда процентное содержание магнетита в нефти определяется по формуле: 

А/ Мрезо4 • 1 00 % = (2,5- 1 0' 5 /92) 1 00 % = 2,7- 1 0' 5 % 

Следует отметить, что полученные значения по порядку величины совпадают с 
определенным в лаборатории атомно-абсорбционного анализа Института органической 
химии УрО РАН (2Т0 ' 5 % Ре). 

Намагниченность (М'р е зо 4 ) одной частицы с массой ш ~ ІО' 14 г. составит 92 ІО' 14 
еши/§. Принимая плотность нефти за 1 г/см 3 определим число N таких частиц в 1 г/см 3 

А • р неф™ / М ' р е3 04 = 2 , 5 Т О' 5 / 92' 1 (У 14 = 2,7- 1 0 +7 1 / см 3 

Итак, в одном кубическом сантиметре мы имеем порядка 30 млн примесей из 
оксидов железа общей массой ~ 2,7- ІО" 7 г. Это, с одной стороны, очень маленькая масса 
ферромагнитных частиц , которую определил только СКВИД-магнитометр, но с дру- 
гой, очень большое число потенциальных зародышей для отложений АСП. Содержание 
ферромагнитных примесей в АСПО резко увеличивается по сравнению с нефтью. Это 
следует из анализа магнитных измерений образцов отложений на внутренней поверх- 
ности НКТ различных месторождений. У АСПО Чернавского месторождения Удмур- 
тии ^180-10‘ 6 , а Тевлинско-Руссинского месторождения Когалыма х~513Т0" 6 . Спек- 
тральный анализ этого АСПО определяет наибольшие проценты для трех элементов 
(кроме органики): диамагнетиков хлора и серы с х = -40- ІО" 6 и -20- ІО" 6 соответственно и 
парамагнетика натрия х = 16- ІО" 6 .Очевидно, что натрий находится в АСПО в соедине- 
нии с хлором, а ЫаСІ уже диамагнитен с х = -0,5 -ІО" 6 . Поэтому положительное значение 
магнитной восприимчивости АСПО целиком определяется парамагнетизмом асфальте- 
нов [10] и ферромагнетизмом примесей. На рис.2, приведены петли гистерезиса образ- 
цов АСПО трех Пермских месторождений, измеренные на стандартном вибромагнито- 
метре. На рис. За. представлены кривые намагничивания АСПО того же Гремихинского 
месторождения, а на рис. 36 - кривые размагничивания, свидетельствующие, что 
АСПО обладает остаточной намагниченностью и коэрцитивной силой. Эти результаты 
получены на высокочувствительном СКВ ИД - магнитометре. 

Проведем анализ этих результатов по аналогии с расчетом, описанном выше. 
Приведенные зависимости демонстрируют, что магнетизм АСПО определяется диа- 
магнитными, парамагнитными и ферромагнитными компонентами. Величина коэффи- 
циента А, определяющего ферромагнетизм АСПО, на три порядка выше чем у нефти. 
Значение удельной намагниченности образцов составляет величину 0,02 ети/§. Для 
оценочных расчетов используем плотность для АСПО тоже 1 г/см 3 . Тогда количество 
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Рис. 2. Петли гиситерезиса трех АСПО с Пермских месторождений: 
1 - Рассвет; 2 - Гожон; 3 - Змеевское 


частиц Ре 3 0 4 с массой порядка 10 ’ 14 г. и намагниченностью одной частицы92 10 -’ 4 
еши/§ в одном кубическом сантиметре составит: 

N = 0,02/92 ІО' 14 = 2.2 • 10 10 1/ см 3 , 

что в тысячи раз выше их концентрации в нефти. 

Отметим важный результат, следующий из рис .2 и рис.З. - ферромагнитная фа- 
за в АСПО имеет коэрцитивную силу и является постоянным магнитом. Можно сказать 
что АСПО-магнитопласт. 

Известно, что частицы ферромагнетика, достаточно малых размеров, являются 
однодоменными (что тоже самое - постоянными магнитами). Еще в 1935 г. Готтшальк 
[11] показал, что коэрцитивная сила природного магнетита резко увеличивается при 
уменьшение размера. Учитывая приведенные данные, можно предположить, что в до- 
оываемои нефти ферромагнитные примеси являются микромагнитами. Тогда один из 
возможных и наиболее эффективных механизмов влияния магнитной обработки нефти 
на процесс отложений АСПО на стенках НКТ можно описать следующим образом. 

Без магнитного индуктора, эти частицы, за счет магнитного притяжения прили- 
пают к поверхности ферромагнитной НКТ вдоль всей ее длины и являются центрами 
кристаллизации АСПО. 
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Рис. 3 . Кривые намагничивания - А и размагничивания - Б АСПО 
Гремихинского месторождения 
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В магнитном индукторе, поставленном в НКТ, на пути потока нефти задержива- 
ется значительная часть этих частиц. Это подтверждается измерениями кривых намаг- 
ничивания отложений уже на МИОНах, бывших в эксплуатации. На рис. 4. приведены 


кривые намагничивания отложений , кривая А - для МИОНа с глубины 900 м и кривая 
В - для МИОНа с глубины 2000м. Удельная намагниченность образцов составляет ве- 


личину порядка 20 ети/§! Спектральный анализ этих отложений показывает наличие 
железа (37% -2000 м. и 30% -900 м.). 

При большом наполнении частиц в индукторе, происходит их срыв потоком 
нефти. В этом случае, ориентированные определенным образом и дополнительно под- 
магниченные полем индуктора, они объединяются в крупные фрагменты, не обладаю- 
щие результирующей намагниченностью, т.к. частицы притягиваются друг к другу раз- 
ноименными магнитными полюсами и образуют замкнутую магнитную систему. Эти 
фрагменты уже не притягиваются к поверхности ферромагнитной НКТ и поднимаются 
вверх в потоке нефти. При этом они являются зародышами кристаллизации асфальте- 
нов, смол и парафинов во всем объеме нефти, а не только на поверхностях оборудова- 

НИ я. 
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Рис. 4. Кривые намагничивания отложений на МИОНах: 
А - на МИОНе с глубины ~ 900 м; 

Б - на МИОНе с глубины ~ 2000 м 
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Очевидно, что эффективность этого процесса определяется величиной напря- 
женности магнитного поля и временем его воздействия на нефть, а это в первую оче- 
редь зависит от ее вязкости. Чем больше вязкость, тем большие поле и время необхо- 
димы для такой кристаллизации в объеме. 

Следует отметить, что и парамагнитные асфальтены в нефти играют роль цен- 
тров кристаллизации парафиновых углеводородов [10]. В магнитном поле они образу- 
ют центры кристаллизации более эффективно, слипаясь за счет магнитного взаимодей- 
ствия в объеме жидкости. 

Но при наличии в продукте ферромагнитных частиц кристаллизация АСП будет 
в основном определяться вышеописанным процессом. 

В литературе описывается еще ряд механизмов воздействия магнитного поля на 
жидкость, содержащую ионы. Рассмотрим некоторые из них. 

Эффект Холла. При наложении магнитного поля на электролит, в котором про- 
текает слабый электрический ток, происходит поперечное разделение ионов разного 
знака, что сопровождается возникновением разности потенциалов на измерительных 
электродах. 

Силы Лоренца. При магнитной обработке возникновение электромагнитных по- 
лей является следствием перемещения воды в магнитном поле. В этих условиях возни- 
кает противоположное циклотронное перемещение катионов и анионов под воздейст- 
вием сил Лоренца, обуславливающих движение ионов вокруг силовых линий магнит- 
ного поля с определенной частотой. 

Эффект Гартмана. При протекании токопроводящей вязкой жидкости между по- 
люсами магнита профиль скоростей потока меняется - параболическая форма искажа- 
ется. 

Все эти эффекты являются по интенсивности воздействия на водяные растворы 
эффектами второго порядка малости и влияние их на скорость отложения АСП не мо- 
жет быть существенным. Сами авторы, использующие в своих работах эти механизмы, 
пишут, например, что «действие сил Лоренца на ионы раствора, перемещающихся в 
магнитном поле, очень невелики» [12]. Оценки показывают, что на длине свободного 
пробега заряженной частицы в жидкости (~10 7 м) силы Кулона превышают силы Ло- 
ренца на четыре пять порядков. 

РЕЗУЛЬТАТЫ ИСПОЛЬЗОВАНИЯ МИОНОВ ПРИ ДОБЫЧЕ НЕФТИ 

В настоящее время в ЦНТ «ЛАНТАН» разработаны малогабаритные высокоградиент- 
ные магнитные индукторы обработки нефти (МИОН) на основе постоянных магнитов 
неодим-железо-бор и самарий-кобальт. 

На рис. 5 схематически показано устройство трубного МИОНа (Патент 
№2198849 от 10.04.2001г.). В ферромагнитной трубе 1 с цилиндрическим пазом разме- 
щены кольцевые магниты 2, каждый из которых образован двумя полукольцами намаг- 
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ниченными радиально противоположно и обращенными друг к другу разноименными 
полюсами. Внутренняя труба 5 для потока жидкости выполняется из немагнитного ма- 
териала. Таким образом устроены трубопроводные МИОНы, имеющие маркировку «Т- 
Ду», где «Ду» обозначает условный диаметр трубопровода («ЦНЪЛантан» производит 
МИОН Т с диаметром от 60 до 200 мм). Такое же устройство имеет и скважинный муф- 
Т °выи ( СМ ) МИОН для НКТ 73. МИОН-СМ73 выполняется в виде традиционной муф- 
ты, соединяющей отдельные трубы, по группе прочности «Е», имеет соответствующие 

резьбы с двух сторон по ГОСТ 633-80 и может устанавливаться в любом месте колонны 
НКТ. 

На рис.6, показано устройство скважинного штангового (СШ) МИОНа (Патент 
№2235690 от 30.10.2002г.). На традиционных штангах 1 по ГОСТ 13877-80 располага- 
ются кольцевые магниты 2, каждый из которых образован двумя полукольцами намаг- 
ниченными радиально противоположно. Магниты запрессовываются в полиамидные 
скребки - центраторы 4, применяемые для очистки НКТ - позиция 5. Сейчас выпуска- 
ется штанговые МИОНы с диаметром штанг 19 (МИОН-СШ 19) и 22 (МИОН-СШ 22), 
которые могут устанавливаться в любое место штанговой колонны. 

На рис. 7. приведена схема установки магнитных индукторов в скважине. Пер- 
вый МИОН-СМ 73 устанавливается на приеме насоса, а второй на границе начала от- 
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Рис. 6. Схема устройства МИОНа штангового типа 



Рис. 7. Принципиальная схема установки магнитных индукторов в скважине 
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ложения АСП. Если используется электроцентробежный насос (ЭЦН), то второй 

МИОН тоже СМ 73, а если штанговый глубинный насос (ШГН) - то второй МИОН СШ 
19 (или 22). 

Первый МИОН был установлен 25.03.2000 года на скважине № 157 Чернавского 
месторождения НГДУ «Удмуртторф». После установки МИОНа горячую обработку 

провели 01.03.2003 - наработка составила 1100 суток. Всего до 12.09.2004 года провели 
6 горячих обработок. 

С июня 2001 года, после опытной эксплуатации, началось внедрение МИОНов в 
ОАО «Удмуртнефть». По данным на 01.01.05 г. на этом предприятии 425 магнитных 
индукторов типа СМ-73 и СШ-19 установлены на 235 скважинах различных месторож- 

дений. Общая наработка составила 133928 суток. За этот период число горячих обрабо- 
ток уменьшилось на 1505. 

В ОАО «Белкамнефть» на 52 скважинах установлено 92 МИОНа СМ-73 и СШ- 
1 9. При общей наработке в 29635 суток на 442 уменьшилось число горячих обработок. 

Результаты применения МИОНов на Тевлинско-Руссинском месторождении Ко- 
галыма приведены в таблице. 

В настоящее время МИОНы эффективно работают еще на ряде нефтедобываю- 
щих предприятий в Удмуртии, Перми, Бузулуке, Ульяновске, Татарстане, Когалыме 
Лангепасе, Казахстане. 
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«Из искры возгорится пламя» 
А.С. Пушкин 

7-9 сентября 2005 г. в г. Берлине (Германия) на базе Федерального института материа- 
ловедения проходила X Европейская конференция по снижению горючести полимеров. 
В этой конференции приняли участие 250 ученых из 23 стран мира (Германии, США, 
Японии, Китая, Англии, Италии, Франции, Израиля, Бельгии, Испании, Тайваня, Гол- 
ландии, Швейцарии, Ирана, Венгрии, Южной Кореи, Австрии, Польши, Швеции, Нор- 
вегии и России), представлявших 1 00 исследовательских центров. 

В состав Оргкомитета конференции входили Ог. ВегпЬагО 8 Нагіе 1 (председатель 
Оргкомитета, Германия), Иг. 8ег ё е ВоигЬщоІ (Франция), Ог.МісЬеІ Бе Вгаз (Франция), 
РгоГ Оюѵапш Сатіпо (Италия), РгоГ Пепе ОеІоЬеІ (Франция), Эг. МГгеу \Ѵ.Оі1шап 
(США), РгоГ А.КісНагё Ноггокз (Англия), РгоГ МепасЬеш Геѵѵіп (США), РгоГ Оуг§у 
Магозі (Венгрия), РгоГ Оеппіз Ргісе (Англия), РгоГ Лапср \Ѵап ё (Китай), РгоГ. СНагІез А. 
\Ѵі1кіе (США), Г.Е.Заиков (Россия) и Мгз. КозхѵііНа Ииззег (ученый секретарь). 

На открытии конференции выступили проф. МапГгеб Неппеске (Президент Фе- 
дерального института материаловедения) и проф. В.ЗсЬагіеІ. Проф. Неппеске сказал, 
что девизом Федерального института материаловедения в его работе является: «Безо- 
пасность и надежность в химических и промышленных технологиях». Участники кон- 
ференции приняли этот лозунг и для своей работы. 

Программа конференции включала 11 секций. На первой секции профессор М. 
Ьеѵѵт выступил с обзорным докладом «Снижение горючести полимеров: прошлое на- 
стоящее и будущее». К.Е.Еуоп (США) рассказал о проблемах термохимии при горении 
твердого тела. 

На второй и третьей секциях («Снижение горючести полимеров») были рас- 
смотрены такие проблемы, как вспучивание полимера при горении (Б.Бияиезпе, 
8.ВигЬі§оІ, КЮеІоЬеІ, Франция), механизмы замедления горючести полибутилентереф- 
талата с помощью красного фосфора (ІГВгаип, В.8Ьаг1е1, Германия); термическая ста- 
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бильность и горение полипропилена с включениями фосфата аммония, микрокапсули- 
рованною в полиуретановой оболочке (З.Оігаиб, А.Сазггоѵіпсі, І.Ѵгошап, О.Сашіпо, 
Италия) Кроме того, были обсуждены вопросы, связанные со снижением горючести 
полимеров с помощью новых методов (М.Иоегіп§, З.ЗеіЬоЫ, И.АгпоШ, М.СіезіеІзкі, 
И.ЗЮггег, ІШіИгісЬ, ВЧизІ, Н.КеІІег, Ѵ.АНзІасК, ІКЖ Запёіег, К. Ргег, В.ЗсЬаПеІ, Гер- 
мания) а также снижение горючести и улучшение механических свойств твердых эпок- 
сикомпозитов с введенными антипиренами (В.К.КапёоІа, Р.МуІег, А.К.Ноггоскз, М. Е1- 
Насіісіі, И. Віаіг, Англия). Использование эпоксидных смол пониженной горючести для 
отливки жидких композитов под давлением (К. Регег, ЖЖ Запсііег, V. АНзГасіі:, 
Т.Нойтап, И.РозріесЬ, М.СіезіеІзкі, М.Иогт§, Германия). Заключительные доклады, 
представленные на этих секциях, были посвящены исследованию термического пове- 
дения антипиренов на основе ковалентносвязанного фосфора (И.Ргісе, Ь.К.СипІіВе, 
КЛ.ВиНеП-РугаЬ, Т.К.НиІІ, О.ГМіІпез, ІК.ЕМоп, В.ЕНшй, РЛозерЬ, Англия) и сниже- 
нию горючести полиакрилонитрилов, ковалентно модифицированных фосфорсодер- 
жащими группами (І.К.ЕЬёоп, Ѵ.Е.Сгоок, В.ГНшй, РЧозерЬ, Р.\Ѵушап, Англия). 

Четвертая секция «Промышленное применение, регулирование и стандартиза- 
ция» включала 4 доклада. В первом докладе была приведена информация о поведении 
различных потребительских товаров при пожаре 5-ой степени сложности( 
М.М.НігзсЫег, США). Кроме того, в докладах была представлена информация об эпок- 
сидных слоистых пластиках пониженной горючести, содержащих антипирены на осно- 
ве фосфора и неорганические наполнители (8.Ѵ. ЬеѵсЫк, М.Висгек, З.ОазЬеѵзку. 
США), о новых экологически чистых антипиренах на основе фосфора (З.НогоШ. 

0. 3сйакег, Германия), а также об использовании редкого типа аморфной двуокиси 
кремния в качестве добавки для снижения горючести термопластиков. 

На 5-ой секции «Текстильные изделия» рассматривались вопросы, связанные с 
заменой галогенсодержащих антипиренов в текстильной промышленности 
(А.К.Ноггоскз, Р.ТОаѵіез, В.К.КапёоІа, Англия), с оптимизацией структуры нанокомпо- 
зитов на основе полиуретан/глина, используемых в качестве антипиренов в тканных 
покрытиях (М.КосЬегу, М.ТакІак, З.ОігаисІ, М.Ьезѵапсіоѵѵзкі, Е.Оеѵаих, З.ВигЬщоГ 
М.АІехапёге, Р.ОіЪоіз, Франция) Кроме того, в докладах, представленных на этой сек- 
ции, обсуждались современные достижения в области снижения горючести текстиль- 
ных изделий (8.ВоигЬі§о1, З.Вияиезпе, Р.РІатЬагб, 8.Ве11ауег, М.КосЬегу, 8.0ігаи<1 

1. Ѵогтап, Франция) в частности, текстильных изделий, применяемых в промышленно- 
сти (Х.ПатЬагб, Б.ЗаіЫ, М.РосЬегу, 8.ВоигЫ§о1, Р.РоиІсЬ, Франция). 

Секция «Новые материалы, новые методы» включала 5 докладов. В первом док- 
ладе обсуждались новые современные методики определения горючести полимеров 
(ГОіІтап, ГЗЫеИз, К.Нагпз ,П\, С.МоШ^ошегу, К. Эаѵіз, США). Большой интерес вы- 
звало рассмотрение таких проблем, как оценка влияния наноглин и тригидрата алюми- 
ния на процессы вспучивания с применением конкалориметра (А.СазІгоѵіпсі, К.Ьуоп. 


182 


ХИМИЧЕСКАЯ ФИЗИКА И МЕЗОСКОПИЯ. Том 7, № 2 


X ЕВРОПЕЙСКАЯ КОНФЕРЕНЦИЯ ПО СНИЖЕНИЮ ГОРЮЧЕСТИ ПОЛИМЕРОВ 


К.\Уакег8, 8.І.8іо1іагоѵ, Т.КазЫ\ѵа§і, О. Сагпіпо, США, Италия), а также разработка но- 
вых антипиренов, не содержащих галогенов, для получения промышленных пластиков 
с помощью оптических методов исследования ХМ.ЕпсЙпег, Германия). Новые подходы 
к методам получения вспенивающихся покрытий и новые методы, способствующие не- 
скапыванию в негорючих термопластиках, обсуждались в докладах, сделанных уче- 
ными из Германии, Израиля и Нидерландов (М.ВагйюІтаі, В.ЗЬагіеІ, Оегшапу; Ѵ.Ваг 
Уаакоѵ, Р.ОеогІейе, К.Вогтз). 

Две секции (7 и 9) обсуждали тему «СО, дым и окружающая среда». На секции 7 
было заслушано 4 доклада. В первом докладе были рассмотрены вопросы оценки сте- 
пени токсичности, возникающей при пожарах различной степени сложности Т.К.НиІІ, 
К.ЬеЬек, К.Т.РаиІ, Н.Ргісе, Англия). Следующие три доклада были посвящены влиянию 
подавителей дымообоазования на горючесть ненасыщенных полиэфирных нанокомпо- 
зитов (8.Ыагаге, В.К.КапбоІа, А.К.Ногтоскз, Англия) и определению токсичности при 
горении электрических кабелей как функции эквивалентности (К.ЬеЬек, Т.К.НиІІ, 
В Рпсе, К.Т.РаиІ, Англия). Кроме того, на различных моделях, включая мебель и ткани 
на основе полиуретанов, рассматривались проблемы влияния кислородного индекса 

при горении (Ь.Тап§е, О.ВизгагсІ, С.ВІипсіеІІ, М.8ітопзоп, Р.Апёегзоп, Н.8ігірр1е, Ни- 
дерланды). 

На секциях 8 и 10 обсуждались нанокомпозиты. В частности, на 8 секции рас- 
сматривались вопросы, связанные с механизмами понижения горючести полимерных 
нанокомпозитов (В. 8сЬайе1, М.ВагйюІтаі, ІГКпоІІ, Германия) с термодеструкцией и 
горючестью полиметилметакрилата с включением микрочастиц А100Н, А1 2 0 3 и нано- 
композитов (А.ЬаасЬасЫ, М.СосЬег, М.РеггіоІ, Е.Ьегоу. ТМХорег-СиезІа, Франция). 

1 акже была представлена интересная и важная информация о кинетическом изучении 
термодеструкции полимерных нанокомпозитов (Г.Е.Заиков, С.М.Ломакин. 

И.Л.Дубникова, С. М. Березина, Россия, и К.КоЛолѵзку, Польша, Оуеоп ё Мап Кіш, 
О.Н.МісЫег, Германия). 

На секции 9 было представлено 5 докладов. В первом докладе обсуждались пра- 
вила и стандартные тесты при пожарах, принятые в Европейском Союзе в настоящее 
время , а также их изменения и перспективы в будущем (в сторону ужесточения) 
(ІН.ТгоігзсЬ). Также была представлена информация о результатах различных тестов 
для проверки безопасности строительных конструкций (Е.АпЮпаШз, Германия). Кроме 
того, обсуждалось моделирование различных сценариев пожаров в жилых зданиях для 
того, чтобы проверить действие новых детекторов дыма и замедлителей горючести 
А.Нойпап. I ермания). Большой интерес вызвало рассмотрение вопроса, касающегося 
применения различных вспененных продуктов с улучшенными параметрами горючести 
(А.ЬогепгеЛі, Италия). 

На секции 10 было представлено 5 докладов. Наибольший интерес вызвала ин- 
формация об образовании нанокомпозитов как компонентов системы для замедления 


ХИМИЧЕСКАЯ ФИЗИКА И МЕЗОСКОПИЯ. Том 7, № 2 


183 


Г.Е.ЗАИКОВ, М.И.АРЦИС, Л.Л.МАДЮСКИНА 


горения (СЬ.АЛѴіІкіе, США), а также информация о наноструктурных замедлителях 
горения (О.Веуег, Бельгия). Часть докладов была посвящена сравнению механизма за- 
медления горения нанокомпозитов в термопластичных и термоустойчивых матрицах 
(А.ТоІсіу, О.Магозі, Р.Аппа, Ы.Тойі, Венгрия). Кроме того, обсуждались получение и 
горючие свойства нанокомпозитов, основанных на эпоксирезинах и органически моди- 
фицированных глинах (А.НаЧ\ѵі§, П.РШг, Германия); физические и горючие свойства 
полимер/кремний и полимер/алюминий нанокомпозитов (Р.Уап§, І.Во§<Запоѵа, 
К.Уп§агб, А.НетЬег§, ОХ.ЫеІзоп, США).. 

Последние 4 доклада были представлены на секции 1 1 . Была рассмотрена ин- 
формация об эмиссии некоторых замедлителей горения в качестве перспективы разра- 
ботки новых замедлителей горения с улучшенными свойствами (8.Кеште1еіп, О.НаЬп, 
ОЧагт, Германия); о правилах, которые приняты в ЕС для урегулирования замедления 
горючести полимеров и реальном соотношении между теоретической и прикладной 
наукой в этой области (О.Визгагё, Англия). Информация о применении трехокиси 
сурьмы в качестве замедлителя горючести для полимерных и композиционных мате- 
риалов (современное состояние и оценка риска) (І-Р.ОаисіеІ, К.Ѵап <3е ѴеЫе, Франция, 
Бельгия) и о замедлителях горения для целлюлозных волокон(К.Ко 2 Іо\ѵзкі, \У1асіука- 
РігЬуІак, Польша) также обсуждалась на этой секции. 

На конференции было также представлено 38 стендовых докладов, посвящен- 
ных различным проблемам горения полимерных материалов и снижению горючести 
полимеров. 

Следующая 1 1 конференция будет проводиться на базе исследовательского ин- 
ститута в г. Болтон (Англия). Организаторами этой конференции являются всемирно 
известные ученые в области горения полимеров профессор Э.Ргісе и профессор 
А.К.Ноггоскз. 
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« Благодарными бывают только пчелы, 
и благодарность их сладкая» 
Шерлок Холмс 
(по замыслу Артура Конан Дойла) 

Как правило, Шерлок Холмс всегда был прав в своих выводах и оценках, применяя де- 
дуктивный метод анализа событий. Но тут его идеи были бы не приемлемы. Наша 
страна, как преемник СССР, Российские ученые, Мировая научная общественность, 
ученики академика Н.М. Эмануэля, его друзья, товарищи, все те, кто его знал и встре- 
чался с ним, сохранили о нем благодарную память. 

Первого октября 2005 года Николаю Марковичу исполнилось бы 90 лет. Россий- 
ская Академия Наук (Р АН), Отделение Химии и Наук о Материалах (ОХНМ), Институт 
Биохимической физики (ИБХФ) им. Н.М. Эмануэля, Химический факультет Москов- 
ского Государственного Университета (МГУ) им. М.В. Ломоносова, Институт Химиче- 
ской физики (ИХФ) им. Н.Н. Семенова и Институт Проблем Химической физики 
(ИПХБ) в Черноголовке (Московская область, Ногинский район) организовали юби- 
лейные мероприятия и конференцию, посвященные 90-летию со дня рождения акаде- 
мика Н.М. Эмануэля. Название конференции («Научные идеи академика Н.М. Эману- 
эля и современная наука») отражает суть событий, т.к. Николаем Марковичем было 
высказано большое количество идей и гипотез, которые, затем, нашли свое прекрасное 
развитие в современной науке. 

Оргкомитет юбилейных мероприятий и конференций возглавил Вице-Президент 
РАН Н.А. Платэ. В него так же вошли академики В.А. Кабанов, С.М. Алдошин, М.В. 
Алфимов, А. И. Коновалов, В.В. Лунин, О.М. Нефедов, М.А. Островский, Л.А. Пирузян 
и А.Е. Шилов, Члены-Корреспонденты РАН К.М. Дюмаев, Г.Б. Манелис, В.Ф. Разу- 
мов, доктора наук С.Д. Варфоломеев, А. А. Попов, Е.Б. Бурлакова, Г.Е. Заиков и к.х.н. 
В.П. Балахнин. 
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Мемориальное заседание Ученого Совета, посвященное юбилею, состоялось в 
МГУ на Химическом факультете 30 сентября 2005 года. С воспоминаниями о Николае 
Марковиче и его роли в развитии Советской науки и воспитании молодого поколения 
выступили академики В. В. Лунин (декан Химфака), Н.А. Платэ, В.А. Кабанов, А.Л. Бу- 
чаченко, И.П. Белецкая, проф. Е.Б. Бурлакова и проф. М.Я. Мельников. 

Следующее заседание состоялось в ИБХФ 3 октября, где и была открыта науч- 
ная конференция. Со вступительным словом на конференции выступил Н.А. Платэ. 
Директор ИБХФ С.Д. Варфоломеев рассказал о жизненном пути и научных достижени- 
ях Н.М. Эмануэля, а д.м.н. М.Р. Лечиницер остановился на значении фундаментальных 
исследований Н.М. Эмануэля для становления клинической онкологии. В обсуждении 
докладов на этом заседании выступили В.А. Кабанов, А.Е. Шилов, А. А. Берлин, К.М. 
Дюмаев, Л.А. Пирузян, Г.Б. Манелис и Л. В. Забелин. 

Второе заседание конференции состоялось в ИПХФ 4 октября. Здесь со вступи- 
тельным словом выступил директор ИПХФ академик С.М. Алдошин. Заседание вклю- 
чало четыре доклада. Н.П. Коновалова посвятила свое выступление антиоксидантам и 
донорам монооксида азота (противоопухолевые эффекты нитроксидов и ГЮ-доноров), 
а доклад Н.А. Саниной и С.М. Алдошина касался информации о новом классе доноров 
N0 (синтез, строение, свойства и перспективы практического применения сера- 
нитрозильных комплексов железа). Свободнорадикальные механизмы возникновения и 
развития вторичного некроза при огнестрельных ранениях были продемонстрированы в 
докладе Г.Н. Богданова, а Л.Д. Смирнов рассказал о фармакологических свойствах и 
перспективных областях клинического применения антиоксидантов гетероароматиче- 
ского ряда. 

В этот же день в ИПХФ состоялось третье заседание, где были рассмотрены 
пять докладов. Большой интерес у слушателей вызвало сообщение Е.Т. Денисова о ра- 
дикальной химии артемизинина. В докладе Л.А. Островской и Д.Б. Кормана обсужда- 
лись вопросы применения нитрозоалкилмочевины в современной противоопухолевой 
химиотерапии. Особенности пероксидного окисления липидов биологических мембран 
в присутствии сверх малых концентраций антиоксидантов обсуждались в сообщении 
Н.П. Пальминой, а доклад Л. Б. Горбачевой был посвящен биохимическим исследова- 
ниям в области химиотерапии. В последнем докладе третьего заседания (Л. К. Обухова) 
были представлены механизмы старения организма и возможности продления жизни. 

Четвертое и пятое заседания состоялись 5 октября в ИБХФ. Семь докладов было 
включено в четвертое заседание. Г.Е. Заиков информировал слушателей о последних 
достижениях в области снижения горючести полимерных материалов и о применении 
нанокомпозитов в качестве антипиренов (вещества, снижающие горючесть полимер- 
ных материалов), а А.А. Попов рассказал о кинетических закономерностях дест- 
рукции напряженных полимеров (реакционная способность напряженных молекул). 
Структурно-динамические особенности межфазных слоев в наполненных полимерах 
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были изложены в сообщении А.Л. Коварского и Т.В. Юшкиной, а доклад ГБ. Париж- 
ского, И.С. Гапоновой, Е.Я. Давыдова и Т.В. Похолок был посвящен изучению меха- 
низмов генерирования стабильных азотсодержащих радикалов в присутствии оксидов 
азота. С большим интересом был заслушан доклад группы авторов (академик Ю.Б. Мо- 
наков, Ю.С. Зимин, В.Д. Комиссаров) о кинетике и механизме озонированного окисле- 
ния кетонов. Последние два доклада этого заседания были посвящены анализу первич- 
ных фотохимических процессов в молекулярных красителях и родственных органиче- 
ских соединениях (В. А. Кузьмин, П.П. Левин, А.С. Татиколов) и исследованиям в об- 
ласти хемилюминесценции (Р.Ф. Васильев). 

Шесть докладов было рассмотрено на последнем (пятом) заседании конферен- 
ции. Доклад Е.Б. Бурлаковой назывался «Биоантиоксиданты - новая страница фарма- 
кологии». В нем говорилось об актуальности данных исследований и о больших прак- 
тических применениях. Последние четыре доклада этого заседания были посвящены 
жидкофазному окислению в микрогетерогенных системах (О.Т. Касаикина), теоретиче- 
ским и прикладным аспектам химии пространственно-затрудненных фенолов (В. Б 
Вольева, И.С. Белостоцкая, Н.Л. Комиссарова и Г.А. Никифоров), принципам стабили- 
зации субклеточных структур функционированием (В.Н. Лузиков), механизмам шапер- 
ноподобной активации (Б.И. Курганов) и новым молекулярным механизмам, объяс- 
няющим старение организмов (А.М. Оловников). 

По всем докладам было проведено широкое обсуждение. Работа конференции 
показала, что многие идеи, высказанные Николаем Марковичем Эмануэлем уже нашли 
свое практическое применение, а часть его мыслей и гипотез успешно развивается в 
исследованиях его учеников, а так же многих ученых, работающих в ключевых иссле- 
довательских центрах различных стран и континентов. 
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ОБ ІСОС’2005 

В период с 1 1 по 14 июля 2005 г. состоялась Пятая Международная Конференция 
по внутрикамерным процессам и горению в установках на твердом топливе и 
ствольных системах (ІСОС’2005). Конференция была посвящена 170-летию со Дня 
Рождения Дмитрия Ивановича Менделеева, внесшего огромный вклад в становление 
современной химической науки, развитие технологии производства порохов и 
взрывчатых веществ, а также техники измерений. 

Целью конференции был научный анализ современного состояния и проблем 
развития исследований в области горения конденсированных систем, технологии их 
получения и утилизации, внутренней баллистики и испытаний технических устройств, 
использующих твердые топлива, пороха, взрывчатые вещества и пиротехнические 
системы в качестве источника газообразования и энергии. 

На участие в конференции подали заявки ученые из Москвы, Санкт-Петербурга, 
Тулы, Ижевска, Казани, Уфы, Перми, Челябинска, Новосибирска, Томска, Бийска: 
сотрудники ведущих проектных и научно-исследовательских институтов, таких как 
Московский институт теплотехники, ФЦДТ «Союз», ФГУП «ФНПЦ «Алтай», ФГУП 
«ГНПП «Сплав», Института Химической Физики РАН им.Н.Н.Семенова, Институт 
проблем химической физики РАН, БГТУ им.Д.Ф.Устинова, МВТУ им.Н.Э.Баумана, 
КГТУ им.А.Н.Туполева и других. На конференции были представлены более 50 
докладов, отражающих основные направления развития научных исследований в 
области горения конденсированных систем и механики твердого тела. В работе 
конференции приняли участие представители Китая, Кореи, Японии. 

Была отмечена необходимость дальнейшего повышения энергетических 
характеристик топлив. В качестве способов достижения этих целей наряду с поиском 
новых энергетически перспективных материалов, используемых при производстве 
топлив, отмечалась необходимость повышения гомогенности смесевых топлив для 
обеспечения стабильных высоких энергетических характеристик топлив. В качестве 
основного пути достижения этого были выделены технологии получения и 
использования кластерных материалов. 

Начиная с этого номера журнала, в нем будут опубликованы материалы докладов 
представленных на конференции. 

Оргкомитет конференции 
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Научно-образовательный центр химической физики и мезоскопии. 
Удмуртский научный центр УрО РАН, Россия, Ижевск 

1. КРАТКИЕ СВЕДЕНИЯ ИЗ БИОГРАФИИ. ПЕРЕЧЕНЬ СДЕЛАННОГО 

Д.И. МЕНДЕЛЕЕВЫМ В РАЗНЫХ ОБЛАСТЯХ 

О гениальном ученом Дмитрии Ивановиче Менделееве написано очень много. Есть 
аналитические обзоры его деятельности. Вместе с тем мало, кто рисковал внедряться в 
его творческую кухню, поскольку это чревато спекуляциями и домыслами о его твор- 
ческом мышлении. Тем не менее, искушение, в связи с последними исследованиями 
работы человеческого мозга, очень велико. Действительно, как рождаются и становятся 
гениями? Этот извечный вопрос гносеологии мучает многие поколения ученых. Можно 
предположить, что гены родителей и окружающая обстановка в детстве формируют 
будущие устремления гениев. А можно ли смоделировать их появление? Как функцио- 
нирует мозг гениев? Какие процессы протекают в их организмах? При моделировании 
обычно учитывают начальные условия. В данном случае, кроме генов родителей, пита- 
ния мамы, нужно учитывать время и место рождения, поскольку влияние поля в мо- 
мент рождения нельзя сбрасывать со счетов. Оно тоже может влиять на зарождающие- 
ся процессы уже в самостоятельном организме. Итак, краткие биографические сведения 
Д.И. Менделеева. 

Дмитрий Иванович родился 27 января (по старому стилю) 1834 года в г. Тоболь- 
ске семнадцатым последним ребенком в семье директора Тобольской гимназии и учи- 
лищ Тобольского округа. Он рано лишился родителей: отец умер, когда ему было 13 
лет, а мать умерла в 1850 году. Большое влияние на развитие его природных способно- 
стей оказала мать, Мария Дмитриевна. Она была умной и очень энергичной женщиной, 
владелицей небольшого стекольного завода и делала все, чтобы дать младшему сыну 
хорошее образование, устроив его сначала в Тобольскую гимназию, а затем в Главный 
Педагогический институт в Петербурге. Вот как пишет Менделеев в посвящении к од- 
ной из своих работ о своей матери: «Это исследование посвящается памяти матери ее 
последышем. Она могла его возрастить только своим трудом, ведя заводское дело; вос- 
питывала примером, исправляла любовью и, чтобы отдать науке, вывезла из Сибири, 
тратя последние средства и силы. Умирая, завещала: избегать латинского самооболь- 
щения, настаивать в труде, а не в словах и терпеливо искать божескую или научную 
правду, ибо понимала, сколь часто диалектика обманывает, сколь многое еще должно 
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узнать, и как при помощи науки, без насилия, любовно, но твердо устраняются пред- 
рассудки и ошибки, а достигаются: охрана добытой истины, свобода дальнейшего раз- 
вития, общее благо и внутреннее благополучие. Заветы матери считает священными Д. 
Менделеев». В Главном Педагогическом институте развитию способностей Менделеева 
к наукам способствовали такие знаменитые учителя, как Остроградский (математика), 
Ленц (физика), Воскресенский (химия), Куторга (минералогия). Возможно воспомина- 
ния юности и лекции Воскресенского и Куторги повлияли на его выбор темы диплом- 
ной работы, основанной на проведенных им анализах ортита и пироксена. Работа назы- 
валась «Изморфизм в связи с другими отношениями кристаллической формы к соста- 
ву». Эту работу впоследствии Менделеев назовет одной из «предтеч» открытия Перио- 
дического закона. Изоморфизму он дал такое определение: «Сходство кристаллических 
форм по причине одинаковости атомного строения и одинаковости атомных атмо- 
сфер». Вуз он закончил с золотой медалью и через год защитил магистерскую диссер- 
тацию по химии «Об удельных объемах». После защиты он становится доцентом Пе- 
тербургского университета, читает лекции по теоретической, а затем органической хи- 
мии. В возрасте 25 лет он был отправлен в двухгодичную командировку в Германию, в 
Гейдельберг, где занимался вопросами капиллярности и поверхностным натяжением 
жидкостей. За годы пребывания заграницей им создан пикнометр для определения 
плотностей жидкостей и открыта абсолютная температура кипения жидкостей (крити- 
ческая температура). В 1861 году после возвращения из Германии он пишет первый в 
России учебник по органической химии, за который получает (по представлению из- 
вестного ученого, химика-органика Зинина) Демидовскую премию. Еще до издания 
учебника из печати выходит его теоретическая статья по органической химии «Опыт 
теории пределов органических соединений», идеи которой впоследствии были развиты 
в других областях химии. В 1865 году Дмитрий Иванович защищает диссертацию «О 
соединениях спирта с водой» на степень доктора химии. Эта работа послужила началом 
гидратной теории растворов Д.И. Менделеева и явилась теоретической основой для 
производства главного национального спиртового продукта - водки «Московской осо- 
бенной», которая была запатентована в 1894 году. Через два года Менделеев становит- 
ся заведующим кафедрой неорганической (общей) химии и начинает работу над одним 
из важнейших трудов своей жизни - учебником «Основы химии», в котором дает обос- 
нование фундаментального закона. Периодического закона. В 1 890 г. Менделеев в знак 
протеста против притеснений студентов покинул университет. С 1892 г. он организует 
Г лавную Палату мер и весов и становится ее управляющим. 

Перечислим еще другие достижения Менделеева. Сюда нужно отнести создание 
гидратной теории растворов, открытие общего уравнения состояния идеального газа 
(уравнение Клапейрона - Менделеева), предложение гипотезы происхождения нефти 
из карбидов тяжелых металлов и принципа дробной перегонки при переработке нефти. 
Он также предложил идею подземной газификации углей, создал российский бездым- 
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ный порох (баллистит). физическую теорию весов и точнейшие методы взвешивания. 
Менделеев оставил след во многих дисциплинах, не только в химии, химической тех- 
нологии и физике. Известны его труды по воздухоплаванию, метеорологии, сельскому 
хозяйству, экономике, народному просвещению. Он интересовался разными сторонами 
жизни Российскоі о общества, но главное, что для него было важным — это наука и при- 
ложение ее в жизнь. Он верил, что «посев научный взойдет для жатвы народной». 

Такая разносторонность и «своеобычность» (слово Менделеева) характерна для 
гениальных людей. 

2. ПРОЯВЛЕНИЕ ГЕНИАЛЬНОСТИ. СВИДЕТЕЛЬСТВА СОВРЕМЕННИКОВ 

Сначала нужно дать определение, кто такой гений. Очень емкое определение приво- 
дится в трактате 16 века известного врача и теолога Хуана Уарте под названием «Ис- 
следование способностей к наукам». По его определению «гений - это человек, кото- 
рый из малых посылок делает верные выводы» . 

Гениальность Менделеева отмечали многие современники, но приведем только 
два свидетельства: химика Чугаева и поэта Блока. 

Ниже приведены выдержки из характеристики Менделеева Л.А. Чугаевым: «Ге- 
ниальный химик, первоклассный физик, плодотворный исследователь в области гидро- 
динамики, метеорологии, геологии, в различных отделах химической технологии и 
других сопредельных с химией и физикой дисциплинах, глубокий знаток химической 
промышленности и промышленности вообще, особенно русской, оригинальный мысли- 
тель в области учения о народном хозяйстве, государственный ум, которому, к сожале- 
нию, не суждено было стать государственным человеком, но который видел и понимал 
задачи и будущность России лучше представителей нашей официальной власти. ... Из 
всех признаков, отличающих гениальность и её проявления, два, кажется, являются 
наиболее показательными — это, во-первых, способность охватывать и объединять 
широкие области знания и, во-вторых, способность к резким скачкам мысли, к неожи- 
данному сближению фактов и понятий, которые для обыкновенного смертного кажутся 
далеко стоящими друг от друга и ничем не связанными, по крайней мере до того мо- 
мента, когда такая связь будет обнаружена и доказана». Александр Блок в своём пись- 
ме дочери ученого — Любови — от 15 мая 1903 г. писал: «Твой папа вот какой: он дав- 
но все знает, что бывает на свете. Во все проник. Не укрывается от него ничего. Его 
знание самое полное. Оно происходит от гениальности, у простых людей такого не бы- 
вает. У него есть всё. Такое впечатление он и производит. При нём вовсе не страшно, 

но всегда не спокойно, это от того, что он всё и давно знает, без рассказов, без намё- 
ков...». 

Однако для подтверждения вышесказанного нужны факты. И самым главным 
подтверждением является создание Периодического закона и предсказание неизвест- 
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ных элементов, которые более, чем на полстолетия опередили строгое обоснование Пе- 
риодического закона с позиций волновой теории атома. Конечно, у Менделеева были 
предшественники, работы которых легли в основание здания «Периодический закон», 
которое он воздвиг. И заслуга Дмитрия Ивановича еще состоит в том, что он сумел за- 
щитить свой приоритет основателя этого фундаментального закона. Лучше самого 
Менделеева о его приоритете не скажешь. Вот его слова: «Утверждение закона воз- 
можно только при помощи вывода из него следствий без него невозможных и неожи- 
даемых и оправдания тех следствий в опытной проверке. Ни де Шанкуртуа, которому 
французы приписывают право на открытие закона, ни Ньюлендс, которого выставляют 
англичане, ни Мейер, которого цитировали иные, как основателя Периодического зако- 
на, не рисковали предугадывать свойства неоткрытых элементов, изменять «принятые» 
атомные веса атомов и вообще считать Периодический закон новым, строго поставлен- 
ным законом природы, могущим охватывать еще доселе необобщенные факты, как сде- 
лано мною с самого начала. . .» 

То, что Д.И. Менделеев сумел защитить свое право первооткрывателя Периодиче- 
ского закона, укрепило его славу гения при его жизни, что нечасто случается. 

Разносторонность и энциклопедичность гениальных людей известна. Однако, если 
проанализировать их деятельность, то выявляется определенная логика их поведения и 
четкая, можно сказать «организованная», направленность их деятельности. Иначе гово- 
ря, есть некая идея, которая гениями апробируется в разных областях жизни человече- 
ства. 

Что же было основным стержнем, определяющим направленность деятельности 
Д.И. Менделеева? Как развивалась его логика мышления? 

3. ОСНОВНОЙ СТЕРЖЕНЬ, ОПРЕДЕЛЯЮЩИЙ НАПРАВЛЕННОСТЬ 

ДЕЯТЕЛЬНОСТИ МЕНДЕЛЕЕВА. ПОПЫТКА ЦЕЛОСТНОГО 

ПРЕДСТАВЛЕНИЯ МИРА 

Можно, конечно, ограничиться высказываниями, которые часто используют в биогра- 
фических очерках о Менделееве, о широте его кругозора, глубине и тщательности его 
научных работ. Однако если выстроить последовательно картину его устремлений и 
исследований, то получается, что Дмитрий Иванович был одержим идеей найти перво- 
основу всего сущего. Основной характеристикой этой первоосновы он выбрал массу и 
заключенную в этой массе энергию. Все его главные работы были непосредственно по- 
священы нахождению этих характеристик и их изменений при превращениях веществ, 
а также классификации элементов и веществ по этим характеристикам. Вместе с тем 
Менделеев живо откликался на потребности науки и промышленности. Поэтому не ос- 
тавляя и не сворачивая с главного пути своего движения в науке и технологии, он пы- 
тается сделать так, чтобы выбранное им направление развивалось в более благоприят- 
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ных экономических, образовательных и политических условиях. Остановимся на неко- 
торых основополагающих работах Дмитрия Ивановича с тем, чтобы проследить эволю- 
цию в его мышлении и логику развития творческого поиска в рамках выбранной им, на 
наш взгляд, идеи. Работа по изоморфизму, по словам Менделеева, вовлекла его «в изу- 
чение более всего химических отношений. Этим она определила многое...» Приведен- 
ное выше определение изоморфизма отражало воззрения ученых того времени. Однако 
Дмитрий Иванович уже в то время говорил, что образование кристаллов зависит от сил 
сцепления между частицами и связывал эти силы с температурой. Вполне логично, что 
в магистерской диссертации он анализирует отношения, существующие между удель- 
ными объемами твердых и жидких тел, «сходственных по составу, форме и свойствам». 
Понятно, что анализируя зависимости между удельными объемами жидкости и того же 
состава твердого тела, появляется желание определить силы сцепления между частица- 
ми в твердом теле при переходе его в жидкость, а затем те же силы в жидкости при пе- 
реходе ее в пар. Именно эту цель он поставил себе во время стажировки в Гейдельбер- 
ге. Там он сделал экспериментальное открытие - существование критической темпера- 
туры, при которой в определенных условиях жидкость моментально превращалась в 
пар. Исследования капиллярности и поверхностного натяжения, проводимые в Герма- 
нии, дали толчок дальнейшим попыткам Менделеева определить силы взаимодействий 
между молекулами, которые он связывал с химическим сродством. Причем химическое 
сродство в то время он рассматривал, как «механическое явление». Основной причиной 
взаимодействия химических частиц Дмитрий Иванович считал явление гравитации и 
при этом ссылался на Ньютона: «Еще Ньютон был убежден, что причина химической 
реакции лежит в простом молекулярном притяжении, обусловливающем сцепление. . .» 

Таким образом, основным стержнем творчества Д.И. Менделеева является его 
стремление целостно представить окружающий мир. И вполне логично в связи с этим 
его дальнейшее исследование в области теории пределов в применении к органическим 
соединениям. Сформулировав эту теорию, Менделеев стал первым теоретиком в облас- 
ти органической химии в России. Согласно этой теории, которая была создана в 1861 
году, все органические вещества подразделяются на достигшие предела, т.е. неспособ- 
ные к реакциям присоединения, и не достигшие предела, стремящиеся при реакциях 
присоединения перейти в предельные состояния. В более широком понимании Менде- 
леев заложил основы реакционной способности химических соединений, определив 
пределы их стабильности. 

В дальнейшем Дмитрий Иванович с легкостью осуществил переход из органиче- 
ской химии, несмотря на очевидные успехи в этой науке, в неорганическую химию. Та- 
кой переход был несложным, поскольку Менделеев всегда или почти всегда проводил 
свои исследования в пограничных областях науки, отчетливо представляя, что только 
иніегрированное знание дает объемное понимание и целостное представление о пред- 
мете изучения. Ступенькой к такому переходу становится упоминавшаяся ранее его 
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докторская диссертация, в которой он заложил основы нового учения о растворах. Ес- 
тественно, занятия неорганической химией и написание учебника «Основы химии» по- 
нудили Дмитрия Ивановича навести определенный порядок в классификации элемен- 
тов. Поэтому работа по написанию этого труда в период 1868-1871 гг. протекала парал- 
лельно с обоснованием и формулированием Периодического закона. Как писал Менде- 
леев: «Эти «Основы» - любимое мое дитя. В них мой образ, мой опыт педагога и мои 
задушевные мысли...» Развивая учение о периодичности, Менделеев, может быть не- 
осознанно, следовал законам диалектики и по сути дела в «Основах химии» проявил 
себя, как «философ в химии» (слова Чугаева). Открытое им явление периодического 
изменения свойств химических элементов в то время не имело физического подтвер- 
ждения. Этим обусловлено, как полагает известный исследователь творчества Д.И. 
Менделеева д.х.н. Б.Н. Трифонов, обращение Менделеева в 1871 году к исследованию 
газов. По-видимому, Дмитрий Иванович полагал, что решение проблемы обоснования 
периодичности находится «в выяснении сил тяготения и изучения свойств среды, их 
передающей» (из работы Д.Н. Трифонова «Д.И. Менделеев. Нетрадиционный подход»). 
Поскольку роль среды играет мировой эфир, а эфир, по мысли Менделеева, это особый 
газ с очень малым весом (ІО' 16 г в м 3 ), то исследования Д.И. Менделеева сводились к 
изучению расширения газов в широком интервале давлений и их сжимаемости. Не- 
смотря на определенное критическое отношение к диалектике и ее законам, Дмитрий 
Иванович в своем творчестве следовал этим законам, в частности, закону отрицания 
отрицания, периодически возвращаясь на новом уровне к тем работам, которые были 
сделаны ранее. Его путь в науке напоминает движение по спирали, так же, как предста- 
вил расположение элементов по атомным весам де Шанкуртуа в виде теллуровой спи- 
рали. Действительно, в 70 — 80 годы 19 века Менделеев вновь возвращается к раство- 
рам. В монографии «Исследование водных растворов по удельному весу», опублико- 
ванной в 1887 году, изложена сущность гидратной теории растворов (в общем понима- 
нии теории сольватации). 

По мнению Менделеева растворы представляют собой «жидкие, непрочные опре- 
деленные химические соединения в состоянии диссоциации». И далее, ... «растворы не 
выделяются в область, чуждую атомистическим представлениям, они входят вместе с 
обычными определенными понятиями в круг тех понятий, которые господствуют ныне 
в учении о влиянии масс, о диссоциации и о газах, в то же время растворы представля- 
ют самый общий случай химического взаимодействия, определяемого сравнительно 
слабыми средствами ...» 

В работе, опубликованной в 1902г. и названной «Попытка химического понима- 
ния мирового эфира», Менделеев возвращается к исследованиям газов: «Уже с 70-х го- 
дов у меня назойливо засел вопрос, да что же такое эфир в химическом смысле? Он 
тесно связан с периодическою системой элементов, ею и возбудился во мне. . . Сперва и 
я полагал, что эфир есть сумма разреженных газов в предельном состоянии. Тепереш- 
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ний мой ответ иной... реального понимания эфира нельзя достичь, игнорируя его хи- 
мизм и не считая его элементарным веществом: элементарные же вещества ныне не- 
мыслимы без подчинения их периодической законности...» Поэтому он ввел в перио- 
дическую таблицу нулевую группу и нулевой период. 

Насколько это правомерно, мы не стали бы сейчас говорить, так как наши позна- 
ния в области пико- и наномира только развиваются. Ясно, что представления Менде- 
леева о мировом эфире в настоящее время не согласуются с современными представле- 
ниями о мире. Вместе с тем некоторые положения из многочисленных работ Дмитрия 
Ивановича, особенно представления о динамичных изменениях веществ при наруше- 
нии условий их формирования, являются основой современной теории самоорганиза- 
ции. 

4. ИДЕИ МЕНДЕЛЕЕВА О ВЛИЯНИИ СРЕДЫ И ЯВЛЕНИЯ 

САМООРГАНИЗАЦИИ 

В общем виде самоорганизация систем объясняется переносом энергии и масс в опре- 
деленном направлении, сопровождаемом чередующимися взаимодействиями, в свою 
очередь, приводящими к упорядочиванию системы. На первый взгляд, новые веяния в 
науке, связанные с наномиром и синергетикой, не имеют никакого отношения к Мен- 
делееву, но это только на первый взгляд. 

Работы Д.И. Менделеева, касающиеся изоморфизма, капиллярности и поверхно- 
стного натяжения, гидратная теория растворов, а самое главное, его убеждение, под- 
крепленное экспериментально, о влиянии среды на большинство процессов, протекаю- 
щих в Природе, предопределили развитие многих направлений в науке, в том числе 
изучение самоорганизующихся (самособирающихся) монослоев и самоассоциации на- 
носистем. При этом форма и размеры образующихся наносистем определяются приро- 
дой (строением и энергетикой) среды и/или носителя, в которой или на котором они 
формируются, а также энергетическими условиями, например, температурой. 

Особый интерес в настоящее время представляют воззрения Менделеева о взаи- 
модействиях молекул воды с молекулами или химическими частицами внесенных в во- 
ду растворяющихся веществ. Последние работы об изменении структуры воды под 
влиянием супермалых количеств высокоактивных веществ (наноструктур) представля- 
ют большую ценность для прогнозирования изменений свойств неорганических кри- 
сталлогидратов и материалов на их основе при введении в них гомеопатических кон- 
центраций активных добавок. Эти воззрения также имеют большое значение для созда- 
ния новых эффективных водорастворимых лекарственных средств. Ниже приведем вы- 
держку из работы, посвященной влиянию на воду сверхмалых концентраций примесей 
(СМК). Эта статья А. Самсонова «Как это может быть?», опубликованная в журнале 
«Экология и жизнь». Поскольку эта работа нам импонирует и соответствует нашим 
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представлениям, мы решили привести из нее следующие выдержки с минимальным из- 
менением текста: 

«Из теории растворов следует необходимость учета взаимодействия растворителя 
и растворяемого вещества. Учет того факта, что внесение примеси в воду может суще- 
ственно изменить объемные и поверхностные характеристики самой воды, способен 
дать многое. Вода всегда имеет примеси. Поэтому вполне возможно, что изменяется 
уникальное сочетание примесных характеристик воды, модифицированной белком или 
иным активным веществом (по примесным компонентам, поверхностному натяжению, 
объемному содержанию микропузырьков и образующих их газов). Тогда, если раство- 
рение примеси вызвало изменение в свойствах воды-растворителя, то при добавлении 
новой воды разводится уже не примесь в воде, а сама вода, в которой первоначально 
проводилось растворение. Смешивается «измененная» и стандартная вода, что приво- 
дит к иному расчету концентрации, и в этом случае, если такого рода расчет работает, 
объяснение эффекта СМК будет вполне классическим. 

Однако при сумасшедшем разведении (10' 3 ° — Ю' !02 ) не останется даже следов 
исходной воды, которой в единице объема было ровно столько, сколько определяется 
числом Авогадро 6,02x1 0" 3 , и потому стоит задуматься о неклассических объяснениях... 

Неклассическое (но не противоречащее теории растворов!) объяснение действия 
эффекта СМК может заключаться в предположении, что примесь вызывает образование 
определенных структур воды во всем объеме или на поверхности, от которой, как пока- 
зывает история развития науки (например, в связи с открытием сверхтекучести гелия), 
весьма нелегко избавиться. Отмеченная в статье В.Е. Жвирблиса страсть воды к ассо- 
циации, структурообразованию сегодня бесспорна. Более того, данные ИК- 
спектроскопии показывают, что существует широкий и почти непрерывный спектр 
структур. Реальная вода представляет набор объемных сетевых структур, в каждую из 
которых входит до сотни молекул воды. . . 

В сильно разбавленных растворах (там, где зарождается эффект СМК) мы можем 
столкнуться с автоволновым процессом, когда структурные ансамбли переходят из од- 
ного выделенного состояния в другое под действием самовозбуждающихся в растворе 
волн неравновесных условий. Периодические реакции Жаботинского— Белоусова сего- 
дня прочно вошли в обиход науки. 

Быть может, в воде запускается некоторая периодическая реакция — циклический 
переход одной структуры воды в другую, в результате которого подобные структуры 
размножаются? Возбудив один или несколько таких процессов в объеме воды, исход- 
ная примесь, видимо, становится уже не нужна. 

Дальнейшее добавление воды в раствор должно приводить к разрастанию процес- 
са самокопирования структур. В результате критическим экспериментом для «струк- 
турного» механизма должно стать изучение зависимости эффекта СМК от температу- 
ры. Возможно, критическим параметром является также время, прошедшее от очеред- 
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ного разбавления до обнаружения эффекта...» Как заметил читатель, автор статьи прак- 
тически полностью зависим от высказанных Менделеевым суждений. 

5. ВОЗЗРЕНИЯ Д.И. МЕНДЕЛЕЕВА И СОВРЕМЕННЫЕ ПРЕДСТАВЛЕНИЯ 
И ГИПОТЕЗЫ 

В настоящее время мы отлично представляем, что периодичность свойств химических 
элементов обусловлена послойным построением электронных оболочек атомов и объ- 
ясняется довольно строго волновой моделью атома. Как известно, волновая модель 
атома и квантовая механика Бора появились, по сути, при "подгонке" условий кванто- 
вания к эмпирическому уравнению Ридберга: 

условие Бора - 2пг = пХ; 

уравнение Ридберга - ѵ/с = Я(1/щ 2 - 1/п 2 2 ), 

где пі и п 2 - числа арифметического ряда, ѵ - частота колебаний, с - скорость света, л. - 
длина волны, Я — постоянная Ридберга. Известно также, что энергия кванта определя- 
ется как в = Ьѵ. 

Если выразить частоту ѵ в уравнении энергии кванта через уравнение Ридберга, то 
окажется, что эта энергия пропорциональна разности обратных квадратов натуральных 
чисел: 


8 = ЬсЯ (1/п, 2 - 1/п 2 2 ) = Ь с Я {(п 2 2 - Пі 2 )/(п, 2 п 2 2 )}. 

Таким образом, энергия кванта, соответствующего определенной линии атомного 
спектра, зависит от произведений натуральных чисел и их квадратов. 

Поскольку уравнения Ридберга и Мозли имеют близкую форму написания, иначе 
функции зависят от квадратов чисел арифметического ряда, являющихся эффективны- 
ми квантовыми числами, которые отражают формы орбиталей в атомах, то можно по- 
лагать, что волновая природа атома непосредственно связана с периодичностью повто- 
рения отдельных числовых значений, характеризующих формы электронных облаков в 
пространстве. Тогда закон отрицания отрицания должен проявиться в таких числовых 
рядах, как квадраты и произведения чисел. Попробуем сопоставить периодическую 
таблицу Менделеева с таблицами квадратов чисел и умножения. Вариантов форм элек- 
тронных облаков четыре, а число их перестановок не превышает 10. Аналогично в пе- 
риодической таблице Менделеева количество периодов и групп не превышает число 10. 
Поэтому принимаем следующие условия для анализа квадратов и произведений чисел: 
квадраты чисел, превышающих 1 0, заменяем суммой цифр, представляющих число та- 
ким образом, чтобы полученное значение числа не превышало десяти. Например, 16 
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представим как 7, 25 - тоже 7, 12 соответствует 3, а 36 - 9. При указанном условии по- 
лучим для квадратов чисел ломаную линию, напоминающую кривую проекции графи- 
ческого отображения уравнения Шредингера. Эта линия многократно повторяется по 
всему полю квадратов чисел и включает всего четыре числа: 1, 4, 7, 9. 

Для произведений, используя тоже условие преобразования чисел, получим набор 
ломаных линий (восемь линий) и одну прямую (при умножении на 9). Представляет 
интерес тот факт, что количество ломаных линий точно совпадает с числом групп в пе- 
риодической таблице Менделеева, а изменение характера этих линий соответствует 
усилению металлических свойств по направлению к первой группе и уменьшение их по 
направлению к восьмой. С первой по четвертую линии видно нарастание чисел, а, на- 
чиная с пятой ломаной линии изменения чисел идут в обратном порядке. Девятая линия 
представляет собой прямую линию с повторением цифры 9. 

Таким образом, анализ произведений чисел в какой-то степени корреспондируется 
с предложенной Менделеевым таблицей, включая нулевую (девятую) группу. Это сви- 
детельствует о единстве Природы, когда периодичность имеет место не только в алфа- 
вите Природы, периодической таблице Менделеева, но и в числовых рядах, построен- 
ных из отдельных повторяющихся кирпичиков - цифр. 

Вполне вероятно, что возмущение среды за счет электронных облаков в атомах 
можно рассматривать, как распространение в пространстве затухающих волн. Однако 
если учесть, что слабые возмущения среды имеют волновой характер с соответствую- 
щей периодичностью, как говорил Менделеев, и не происходит мгновенной диссипа- 
ции соответствующих искажений плотности среды (в фоне), то взаимодействие частиц 
в глубоком вакууме, вероятность столкновений которых практически равна нулю, 
можно объяснить за счет интерференции волн и роста интерференционной энергии, со- 
гласно воззрениям Рюденберга. Объяснение сближения химических частиц и образова- 
ние химических связей возможно расширить в общефилософском аспекте на картину 
образования элементарных частиц за счет интерференции волн в пространстве, пони- 
мая под этим, что поле и излучение имеют определенные флуктуации плотности. 

В связи с этим уместно вспомнить серию работ Я.Б. Зельдовича с соавторами о 
рождении элементарных частиц при взаимодействии гравитационных полей. Естест- 
венно, возникает вопрос о сохранении условий квантования в атоме и в полях (излуче- 
ниях) различной природы, а также о взаимном влиянии дискретных полей. 

Любые тела, вещества и химические частицы имеют свои гравитационные, маг- 
нитные и электрические поля. При этом происходит искажение геометрии и энергии 
пространства в месте нахождения соответствующих тел или частиц и в ближайшем их 
окружении. Созданные флуктуации плотностей в пространстве приводят к изменению 
характера и скорости протекающих между химическими частицами, веществами и те- 
лами процессов. Тогда временные промежутки, устанавливаемые в соответствии с про- 
текающими процессами, будут зависимы от искажений геометрии и энергий простран- 
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ства в окрестности частиц. Это, в свою очередь, определяется характером и энергией 
интерференции распространяющихся «в фоне» (или в мировом эфире, согласно Менде- 
лееву) флуктуаций (волн), порождённых колебаниями частиц, дискретностью потоков 
излучений (полей). 

Интенсивность интерференции определяется квадратами амплитуд колебаний со- 
ответствующих составляющих рассматриваемых объектов. Если придерживаться той 
же общей логики, то образование более сложных систем по сравнению с атомами, 
должно быть обусловлено теми же причинами, что и образование элементарных частиц 
и атомов с учётом соответствующего числа перестановок и с подчинением тем же пра- 
вилам квантования, т.е. энергетические параметры более сложных систем вполне ре- 
зонно связывать с квадратами и произведениями чисел арифметического ряда. Переход 
от микро- к макромиру, от элементарных частиц к агрегатам молекул , потокам роёв 
частиц, а также от процессов на атомном и молекулярном уровнях к процессам в ре- 
альных аппаратах и системах при предлагаемом подходе в общем философском плане 
возможен только при «ступенчатом движении» или системном подходе с применением 
принципов и методов согласования. 

Здесь нужно определить, что процессы самоорганизации обусловлены и большей 
степени влиянием среды и периодических возмущений, иногда сравнительно малых, но 
вполне достаточных для установления определённой периодичности чередования 
структурных единиц в системе. 

Если рассматривать образование химической связи как результат интерференции 
волн - электронных облаков, понимая под этим, что рост электронной плотности в 
межъядерном пространстве взаимодействующих атомов сопровождается увеличением 
интенсивности колебаний и «рождением нового осциллятора», то нужно признать из- 
менение возмущений пространства, которое, в свою очередь, приведёт к направленно- 
му искажению геометрии системы в данной области пространства. Изменением вклада 
интерференционной энергии в образование химических связей можно объяснить раз- 
ницу в их природе. Увеличение вклада интерференционной энергии при взаимодейст- 
вии химических частиц может привести к разделению зарядов и образованию ионов, а 
также к изменению характера взаимодействия - росту зарядового и заряд - дипольного 
взаимодействии. Это в ряде случаев определяет перераспределение, а иногда перенос, 
зарядовых плотностей, с изменением структуры системы. 

Возникновение за счет роста энергии интерференции в определенном направле- 
нии перемещений заряда с образованием определенного вида химических связей спо- 
собствует образованию цепочечно-пачечных, плоскостных, цилиндрических (тубуляр- 
ных) и клубковых (глобулярных) структур в зависимости от количества и природы пер- 
воначально образовавшихся фрагментов ("зародышей"). При анализе процессов ради- 
кальной сополимеризации и полимеризации, например, показано, что геометрия расту- 
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щих цепей уже определена на пятой стадии процесса, когда восемь атомов углерода 
соединены ковалентными связями. 

Значительную роль при рассмотрении протекания направленных процессов полу- 
чения и превращений наноструктур, наночастиц и наносистем, процессов, связанных с 
самоорганизацией, играют структурированные и самоорганизующиеся среды. Здесь 
мысли Менделеева о влиянии среды на силы сцепления молекул и атомов, его воззре- 
ния о возможных искажениях в структурах частиц и их свойствах при воздействии на 
них мирового эфира начинают приобретать новое значение. 

В заключении следует сказать, что творческое наследие гениального ученого- 
энциклопедиста Дмитрия Ивановича Менделеева еще требует кропотливого изучения, 
особенно, в приложении к новым реалиям в науке и в промышленности для того, что- 
бы, как говорил Менделеев, "посев научный взошел для жатвы народной". 
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АННОТАЦИЯ. Осуществляется проверка гипотезы влияния неравномерного рас- 
пределения по объему основных компонентов баллиститов, которые обычно счита- 
ются гомогенными композициями, на реализацию локально-нестационарного меха- 
низма горения. Произведен расчет межмолекулярного взаимодействия нитратов, 
рассмотрена решеточная модель истинного раствора мономер-полимер и модель ка- 
пиллярно-пористого тела. 

Сделан вывод об отсутствии термодинамических факторов, препятствующих 
взаимному неограниченному растворению компонентов. Решена задача о проникно- 
вении растворителя в индивидуальную частицу НЦ. Сделан вывод: неравномерное 
распределение компонентов является задающей функцией (задающим полем) в реа- 
лизации локально-нестационарного горения баллиститного топлива. 


ВВЕДЕНИЕ 


Впервые сведения о реализации локально-колебательного механизма горения балли- 
ститных порохов и смесевых топлив были изложены Г.Г.Шелухиным в 1965 г. на кон- 
ференции в Институте химической физике, что вызвало оживленную дискуссию [1]. 
Позднее К. И. Синаев предложил термин «очагово-пульсирующий механизм» приме- 
нительно к горению пиротехнических составов, а затем и нитроглицериновых порохов 
[2]. Для объяснения этого феномена рассматриваются две гипотезы: воздействие ло- 
кальной неоднородности компонентного состава топливной массы и проявление ло- 
кальной неустойчивости процесса горения. В пользу второй гипотезы в чистом виде 
могли бы свидетельствовать данные о горении гомогенных веществ. К таковым часто 
относят базовые баллиститные композиции нитроглицерин - нитроцеллюлоза (НГЦ- 
НЦ). Однако такое представление требует обсуждения. 

Известны экспериментальные данные о влиянии на скорость горения балли- 
ститов сорбционной способности исходной НЦ и существование связи между измене- 
нием скорости горения от объемного соотношения полимер-растворитель и соответст- 
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вующим изменением давлений паров растворителя над системой. Таким образом, ско- 
рость горения оказывается зависящей от термодинамического сродства растворителя 
(НГЦ) к полимеру (НЦ). 

Существовали два различных представления о фазовом составе системы НГЦ- 
НЦ. По одному из них (А.С. Бакаев и сотр.) система не является однофазной и может 
находиться в устойчивом состоянии до массовой концентрации растворителя 0.5, рас- 
сматриваемой как предел термодинамической устойчивости системы. По второму 
(И.В.Тишунин) пластифицированная масса НЦ - истинный раствор полимера. Оконча- 
тельное суждение осложнялось тем, что теоретически возможен вариант системы в ви- 
де наполненного жидкостью капиллярно-пористого тела, и изотерма сорбции может 
совпасть с зависимостью, которая трактуется как зависимость упругости пара от кон- 
центрации компонентов истинного раствора. 

Требовалось исследование термодинамической совместимости компонентов 
системы НГЦ-НЦ, изучение кинетики проникновения НГЦ в нитроцеллюлозу и выяс- 
нение вопроса о степени дисперсности поглощенного НГЦ, необходимого для установ- 
ления состава системы. Тогда становится возможным определение распределений ком- 
понентов по объему топлива и его влияние на локальные характеристики процесса го- 
рения. Актуальность проблемы сохраняется [3] и ниже освещаются эти вопросы. 

О МЕЖМОЛЕКУЛЯРНОЙ СВЯЗИ КОМПОНЕНТОВ БАЛЛИСТИТНОГО 
ТОПЛИВА 

Вывод термодинамических функций для пластифицированных нитроглицерином масс 
нитроцеллюлозы в терминах межмолекулярного взаимодействия представляет извест- 
ные трудности. Из экспериментов известно, что практически вся энергия взаимодейст- 
вия молекул приходится на связь нитратных групп, причем ее величина намного пре- 
восходит среднюю энергию ван-дер-ваальсовых сил. По температурной зависимости 
давления паров над чистым растворителем и над пластифицированной НЦ можно оце- 
нить энергию, приходящуюся на одну связь, в 30 - 40 кДж/моль. Это свидетельствует о 
высокой степени ассоциации растворителя и прочности связей при сольватации. 

Индивидуальную молекулу низкомолекулярных компонентов баллиститного 
топлива, НГЦ и ДНТ и звено НЦ следует рассматривать как совокупность нескольких 
силовых центров, каждый из которых является относительно самостоятельным осцил- 
лятором, образующим мгновенные виртуальные мультиполи. 

Рассмотрим на примере молекулы НГЦ, находящейся в окружении себе подоб- 
ных молекул, пространственное изменение потенциальной энергии взаимодействия. 

На рис.1 изображена стереохимическая структура молекулы НГЦ, построенная 
по данным монографии [4]. Условно два атома С и две нитратные группы совмещены с 
плоскостью чертежа поворотами относительно линии связей, с сохранением валентных 
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Рис.1 . Стереохимическая структура молекулы НГЦ (2 атома и 2 нитратные груп- 
пы совмещены с плоскостью горения, атомы Н не показаны)- 
^СЖО (а) = 125°, ^(ЖО ф) = = 1 1 7° 30', остальные углы близки к тет- 
раэдрическому, составляющему 109°28' 

углов. Атомы Н ради ясности чертежа не показаны. Для представления об объеме мо- 
лекулы проведены сферы со средними ван-дер-ваальсовыми радиусами для каждого 
атома. Штриховой линией изображена сфера со средним радиусом молекулы, который 
получен из предположения о тетраэдрической, упаковке молекул НГЦ. При таком рас- 
чете радиус этой сферы минимален, и тем более легко видеть, что каждая нитратная 
группа может попасть в объем, занимаемый соседней молекулой. Близкое расположе- 
ние активных групп соседних молекул - свидетельство сильной межмолекулярной свя- 
зи. Из рис.1 также следует, что ориентация в пространстве каждой нитратной группы 
обособленной молекулы практически произвольна. 

Возможные варианты межмолекулярной связи: нитратная группа одной моле- 
кулы с нитратной группой молекулы, нитратная группа одной молекулы с ненитратным 
«контактным» участком другой молекулы и, наконец, взаимодействие двух ненитрат- 
ных контактных участков. Такое деление отнюдь не отрицает наличия потенциальной 
энергии связи каждого контактного участка, т. е. отдельной характерной группы атомов 

молекулы со всеми остальными группами атомов второй молекулы (учитывается в рас- 
четах ниже). 

В силу высоких значений экспериментально наблюдаемой энергии связи нит- 
ратных групп, сравнимых по величине с энергией водородных связей, целесообразно 
использовать электростатическую теорию межмолекулярного взаимодействия. Прием- 
лемость модели существенно зависит от точности определения орбитальных углов и 
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распределения зарядов, поэтому приводимый ниже расчет не претендует на большую 
точность, а скорее служит иллюстрацией механизма взаимодействия нитратных групп. 

Будем строить электростатическую модель группы 01чЮ 2 . Воспользуемся 
адиабатическим приближением квантовой механики в определении взаимодействия 
частиц, сильно различающиеся по массе. Колебательным движением ядер пренебрега- 
ем, считаем их закрепленными в центрах колебаний. Связи жесткие. Тогда основной 
задачей станет отыскание центров тяжести полей вероятностей электронной плотности 
с учетом валентного взаимодействия между атомами. 

Трехэлектронное взаимодействие в полуторной связи О N легко учитывает- 
ся в приближенном методе расчета - методе молекулярных орбит [4]. Воспользуемся 
его простейшей модификацией - методом ЛКАО (линейной комбинации атомных ор- 
бит) - в соответствии с точностью расчета. В этом случае волновые функции электро- 
нов, находящихся на молекулярных орбитах, составлены как линейные комбинации 
атомных волновых функций. Согласно принципу Паули на орбите с минимальной 
энергией находятся 2р х электрон азота и 2р х электрон кислорода, на орбите следующе- 
го энергетического уровня - 2р г электрон кислорода. 

С использованием волновых функций многоэлектронного водородоподобного 
атома волновые функции двух электронных орбит с нижними значениями энергии при- 
обретут вид: 


Ѵі = 


0.56 


8л/л 


С \ 


Ѵ а о У 


г м е" г ^ 2а ° зіпѲ ы е' <р * + 


0.83 


8 Ѵл 


/ л 


Ѵ а О у 


г е - 2о '-° 2а о 5 і п ѳ п е ,<Ро 


( 1 ) 


У2 = 


8л/я 


( „ л 


V а о У 


г 0 е 


2 о г о2а 0 


зіпѲ п е 


'Фо 


где 2 - порядковый номер атомного ядра: г, Ѳ, ср - сферические координаты; а 0 - первый 
боровский радиус атома водорода; 0,Ы - индексы, указывающие на принадлежность к 
соответствующим ядрам О и N. 

Для одинарной связи О - N аналогичным путем получена іу, тождественная іу, 
соотношения (1). 

Теперь, зная, что квадрат волновой функции равен плотности вероятности пре- 
бывания электрона в данной точке пространства, нетрудно найти положение центров 
тяжести электронных облаков. Задача эквивалентна отысканию центра тяжести неогра- 
ниченного тела переменной плотности. 

На рис. 2 в аксонометрии показано распределение зарядов в электростатиче- 
ской модели нитратной группы. Было принято, что центры тяжести 1з электронов сов- 
падают с равновесными положениями ядер, а каждое из 2р электронных облаков разби- 
то на две симметричные относительно ядра части со своими центрами тяжести. 
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Рис. 2. Электростатическая модель нитратной группы 

Наконец, электронные облака валентных электронов, находящихся на молекулярных 
орбитах, разбивались на три участка по линии валентной связи координатами ядер. 

Электростатическое взаимодействие входит составной частью в общее выра- 
жение энергии взаимодействия двух нитратных групп, которое, в пренебрежении ма- 
лым индукционным эффектом, имеет вид: 

и _ у д.д, а і а 2 1 В 

Ь іу 1/І,+1/І 2 + ’ (2) 

где ц, , ^ - электрические заряды двух взаимодействующих нитратных групп; г,, - рас- 
стояние между соответствующими зарядами; ос , а 2 - средние поляризуемости 1 и 2-й 
нитратных групп; I) ,І 2 - потенциалы ионизации 1 и 2-й нитратных групп; К. - расстоя- 
ние между центрами тяжести нитратных групп; В - эмпирический коэффициент. Пер- 
вый член выражения (2.10) учитывает кулоновскую энергию взаимодействия, второй - 
дисперсионную, третий - отталкивательную. 

Установлено наличие двух потенциальных ям (рис. 3), различающихся величи- 
ной и, причем меньшей потенциальной энергией обладает вариант, когда плоскости 
треугольников ОЖ) взаимодействующих нитратных групп параллельны (тип I), и не- 
сколько большей, когда они перпендикулярны (тип II) (рис. 4). На расстояниях между 
центрами тяжести групп порядка 3.5 А и выше тепловое движение нарушает взаимную 
ориентацию, которая переходит в случайную. 
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О N О 



Рис. 3. Потенциальная энергия 
взаимодействия двух нитратны 
групп: 1 — по типу I; 2 - по тип 
П; 3 — при взаимодействии с осто 
вом другой молекулы 


Рис. 4. Пространственное располо 
жение потенциальных ям во взаи 
модействии нитратных групп 


Величина В в формуле (2) найдена из предположения о том, что низкомолекулярные 
нитраты в конденсированном состоянии взаимодействуют своими нитратными группа- 
ми по типу I. При этом учитывалось, что к определенной согласно (2) величине и необ- 
ходимо прибавить потенциальную энергию, вызванную остальными атомами двух со- 
седних молекул, согласно Лондону [5]: 


и 0 - 


К 6 


2ц, ц 2 2 
2кТ ^ 2 


+ В2«| 


+ — а.а, - 
2 1 2 І 


М2 

I + ^2 


'бр_і' 

и ) 


( 3 ) 


где ц. 1,2 - суммарные дипольные моменты, отнесенные к центрам тяжести соседних мо- 
лекул; а ід - средние поляризуемости соседних молекул, І 12 - их потенциалы иониза- 
ции; р = 0.28 А. 

Вычисленная при ц = 2.9 дб, К. = 6.5 К., = 2.1 А, I = 920 кДж/моль, а = 22 А 3 и Т 
= 300К, согласно формуле (2.11), и 0 » - 4.2 кДж/моль. Тогда иі * - 35 кДж/моль, В = 

3.6 • ІО 6 кДж А 12 /моль, К] = 2.1 А, и, , * - 25кДж/моль, К„ = 2.6 А (I и II - индексы ука- 
занных двух типов взаимодействия). 

Вычисленная при ц = 2.9 дБ, К. = 6.5 К, = 2.1 А, I = 920 кДж/моль, а - 22 А 3 и Т 
= 300К, согласно формуле (2.1 1), и 0 * - 4.2 кДж/моль. Тогда и, * - 35 кДж/моль, В = 

3.6- ІО 6 кДж А 12 /моль, К, = 2.1 А, и„ * - 25кДж/моль, К„ = 2.6 А (I и II - индексы ука- 
занных двух типов взаимодействия). 
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На рис.4 показаны изопотенциальные кривые, нанесенные в пространстве око- 

ло одной из нитратных групп. Из их рассмотрения следует, что для образования связи 
по типу I необходимо, чтобы одна из соседних молекул обладала свободой вращения 
относительно своего центра тяжести и достаточно свободным вращением нитратных 
групп. По-видимому, все низкомолекулярные нитраты обладают этим свойством. В тех 
случаях, когда обе молекулы имеют заторможенное вращение, реализация взаимодей- 
ствия по типу I затруднена и связь II более естественна. Таким свойством могут обла- 
дать, например, звенья жестких макромолекул НІД. 

Предполагая, что в стеклообразной НЦ элементарная ячейка близка по своему 
строению элементарной ячейке кристаллита целлюлозы по Майеру и Мишу (поправка 
вводится лишь на изменение объема при переходе от Ц к НЦ в соответствии со значе- 
ниями плотностей) и, учитывая устойчивую конформацию а - Б - глюкозных остатков, 
получаем пространственное расположение звеньев НЦ согласно рис. 5. Очевидно, мо- 
жет осуществиться только связь с перпендикулярным расположением групп ОЫСЬ (тип 

И). 

В пластифицированных массах НЦ молекулы растворителя - нитрата могут об- 
разовать более прочные связи со звеньями НЦ (тип I). Разность значений энергии ям I и 

II высвобождается при растворении и служит причиной значительных тепловых эффек- 
тов растворения НЦ. 



Рис. 5. Реализация нитратных связей в НЦ 
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Аналогичные расчеты для взаимодействия групп N02 и 0И02 показывают, что 

величина потенциальной ямы I такая же, как и в случае нитратных групп, а ямы II 

меньше (штриховая линия на рис. 3). 

ТЕРМОДИНАМИЧЕСКИЕ СООТНОШЕНИЯ ДЛЯ ПЛАСТИФИЦИРОВАННЫХ 
МАСС НИТРОЦЕЛЛЮЛОЗЫ 

Воспользуемся обычным приемом, применяемым для установления границ термодина- 
мической совместимости компонентов. Предположим, что осуществляется молекуляр- 
ное взаимное диспергирование НЦ и НГЦ, т.е. образование раствора при всех концен- 
трациях компонентов. Концентрация компонента, начиная с которой химический по- 
тенциал превысит таковой для чистого вещества, будет служить границей совместимо- 
сти. 

Благодаря специфике межмолекулярного взаимодействия нитратных групп, 
характеризуемой резкой анизотропией молекулярного силового поля, к теоретическому 
определению термодинамических функций истинных растворов нитратов необходимо 
привлечь теорию, в которой возможен учет ориентационного эффекта во взаимодейст- 
вии [6]. 

Допущения, лежащие в основе решеточной модели раствора: предполагается, 
что раствор и чистые компоненты обладают квазикристаллической структурой с оди- 
наковым координационным числом; производится учет взаимодействия лишь между 
молекулами, расположенными в непосредственной близости друг от друга, так что 

среднее состояние каждой молекулы является функцией состояния только ближайших 
соседей. 

Потенциальная энергия раствора разбивается на две не зависящие друг от дру- 
га части: конфигурационную, связанную с положениями центров колебаний молекул, 
совпадающих с узлами квазикристаллической решетки, и акустическую энергию коле- 
баний центров масс около узлов решетки. В предположении о зависимости акустиче- 
ской энергии только от числа молекул в растворе нетрудно показать, что связанная с 
ней часть интеграла состояний при постоянном числе молекул постоянна, и все описа- 
ния термодинамических свойств системы приходятся на конфигурационную часть ин- 
теграла состояний. 

Размеры молекул компонентов раствора близки по величине и форме, так что 
осуществляется свободное взаимное замещение ими друг друга в узлах квазикристал- 
лической решетки. Для растворов полимеров, поскольку в качестве термодинамической 
единицы обычно фигурирует сегмент, состоящий из нескольких звеньев полимерной 
цепи, это требование заменяется на условие взаимного замещения сегмента и соответ- 
ствующей ему группы из определенного числа молекул низкомолекулярного раствори- 
теля. 
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В наиболее общем виде учет ориентационного эффекта, свойственного моле- 
кулам нитратов, разработан в теории Баркера [6] применительно к низкомолекулярным 
жидкостям. Ее развитие на полимеры и их растворы приводит к определению термоди- 
намических параметров пластифицированных масс НЦ. Баркер дает зависимости для 
расчета термодинамических функций раствора в «квазихимическом» приближении с 
учетом вероятностей межмолекулярных связей. 


Ф * Р = - кТ ІпО + кТ^Ы^Іп 


А, и 


X 


( ат Х, А ) 


Н«Е = -2ЫкТ^ / Ы А 0 А 1п 




х' 




(А, В = 1,2, ...М) 


(4) 


где Ф, Р, Н, Е - конфигурационные части термодинамического потенциала, свободной 

энергии, энтальпии и внутренней энергии системы, состоящей из N молекул; А, В 

компоненты из числа М компонентов системы; ц, ѵ - виды контактных участков моле- 
кул компонентов; С> А - число контактных участков молекулы компонента А р-го клас- 
са; Б1 а - число молекул компонента А; і! включает в себя члены с р А = ѵ в ; N - полное 
число частиц в системе, О - число конфигураций раствора. Величина О зависит от чис- 
ла конфигураций в предположении атермичности раствора О 0 и энергии взаимодейст- 
вия контактных участков молекул между собой. В свою очередь, О 0 определяется чис- 
лом способов распределения молекул полимера различной длины по узлам квазикри- 
сталлической решетки [6]. 

Представим в качестве модели рыхлого полимера квазикристаллическую ре- 
шетку, часть узлов которой занята звеньями полимера, а часть свободна. Вакантные 
места, «дырки», по терминологии жидкого состояния, могут заниматься отдельными 
звеньями, освобождающими свои узлы, которые становятся вакантными. Так происхо- 
дит обмен местами, аналогичный обмену растворителя с отрезками цепи полимера. 
Учитывая, что кинетической единицей полимера, способной к обмену; является сег- 
мент, состоящий из нескольких звеньев, введем предположение о соединении вакант- 
ных узлов в группы, конгруэнтные сегменту полимера. Таким образом, рыхло упако- 
ванный полимер с термодинамической точки зрения формально подобен раствору, од- 
ним из компонентов которого является плотно упакованный полимер, другим - пусто- 
ты. Ввод низкомолекулярного растворителя в рыхлый стеклообразный полимер приво- 
дит к тому, что часть молекул растворителя стремится заполнить вакантные места, а 
часть раствориться в областях с плотной упаковкой, увеличивая число мест квазикри- 
сталлической решетки. На первый взгляд казалось бы, вначале должно произойти за- 
полнение вакантных мест, а затем растворение. Однако это не так и, если отвлечься от 
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кинетики проникновения растворителя в полимер, то сравнением термодинамической 
вероятности того и другого пути можно указать направление суммарного процесса 
проникновения растворителя в полимер. 

В исходном полимере Н молекул г-мера заполняют гК узлов квазикристалли- 
ческой решетки, т 0 гЫ г узлов остаются вакантными (ш 0 <1). При вводе растворителя в 
полимер число вакантных мест уменьшается, и текущее значение составляет тгК Од- 
новременно увеличивается число мест в квазикристаллической решетке N. 

Очередная молекула растворителя, проникая в полимер, может либо занять од- 


но из вакантных мест, и тогда число конфигураций увеличится в 


шгЫ г 
К +1 


раз в предпо- 


ложении г * я, либо создать новый узел, и тогда число конфигураций увеличится 

N. +1 

раз. 


N. +1 


Вводя относительные вероятности и перейдя от дискретных величин к непре- 
рывным, что возможно при больших числах, получаем следующую систему уравнений: 


сіН= — 1 ^ +1 = сШ, 
N. +1 + т 


, — т 

- бт = =г <йчІ| 

N.. +1 + гп 


( 5 ) 


N$ = N1 + гК + т 
где ш = тгК . 

Начальные условия К = 0; т = т 0 ; К = Ко = гК + т 0 . 

Вводя объемную долю полимера <р = , имеем связь т с объемной до- 

лей полимера в растворе: 

тф= ~ГѴ? +т ° (т ° +1)рі (о 

Очевидно, такой же результат получится, если оперировать в качестве термо- 
динамических единиц сегментами полимера. Тогда число конфигураций изменяет не 
одна молекула растворителя, а группа из а молекул, занимающая дырку, конгруэнтную 
сегменту полимера (в этом случае все значения N нужно разделить на ст). 
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Были приняты следующие данные для расчета: - 39 кДж/моль - энергия взаи- 
модействия нитратных групп по типу I (низкомолекулярные нитраты в связи между со- 
бой и с НЦ); - 4 кДж/моль - энергия взаимодействия прочих контактных участков ме- 
жду собой; - 4 кДж/моль - энергия взаимодействия между нитратными и прочими кон- 
тактными участками между собой; - 26 кДж/моль - энергия взаимодействия нитратных 
групп по типу II (НЦ); дырки имеют, естественно, нулевые значения взаимодействия со 
своим окружением; 2 = 10; ст = 2; г = 1000. 

Результаты расчета для чистых нитратов даны в таблице. 


Нитрат 

Н, кДж/моль 

А8, Дж/моль К 

Мономер с тремя нитратными группами 

-73.5 

-25.4 

Мономер с двумя нитратными группами 

-55.5 

-21 

Мономер с одной нитратной группой 

-38 

- 13.3 

Нитроцеллюлоза ( 1 2% Ы) 

-44’ 

- 18.8* 


* В расчете на моль звеньев. ДЗ - избыточная энтропия, вязанная с ориентацией 
молекул; ДЗ = к (ІпО - 1пО 0 ). 

Аналогичные расчеты для смесей низкомолекулярных нитратов дают энтро- 
пию растворов, практически не отличающуюся от энтропии идеальных растворов, и 
нулевой тепловой эффект при всех концентрациях, т.е. указанные растворы подчиня- 
ются закономерностям идеальных растворов. Это полностью отвечает известным экс- 
периментальным данным по смесям растворителей для НЦ: НГЦ, ДНТ, НГЛ, ДЭГДН 
ТЭГДН и др. 

Выше была отмечена стерическая невозможность у нитратных групп НЦ всту- 
пить в наиболее выгодное взаимодействие вследствие жесткости полимерной цепи. По- 
тенциальная энергия связи имела лишь относительный минимум. Низкомолекулярный 
растворитель НГЦ, обладая подвижностью и свободой вращения около центра тяжести 
молекулы, не встречает препятствий для образования связи как между нитратными 
группами себе подобных молекул, так и во взаимодействии с нитратными группами НЦ 
с абсолютным минимумом потенциальной энергии (тип I). Поскольку движение отно- 
сительно, достаточно условия свободы движения одной молекулы, участвующей в пер- 
вой связи. С переходом на низший энергетический уровень выделяется часть энергии, 
равная половине разности между уровнями, что составляет - 6.5 кДж/моль. 

Благодаря этому термодинамическая вероятность, образования связи нитрат- 
ных групп НЦ и НГЦ в ехр(6500/ЯТ) раз больше вероятности образования связи одно- 
именных молекул. При Т = 300К - в 13.6 раза. Таким образом, механизм сольватации 
НЦ заключается в экранировании одной нитратной группой растворителя одной нит- 
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ратной группы НЦ. Такой взгляд на механизм сольватации был высказан 
Д . И . Гальпериным на основе экспериментальных данных. Здесь получено его термоди- 
намическое подтверждение. 

На рис. 6 показано сравнение тепловых эффектов растворения НЦ в НГЦ по- 
лученных теоретически и экспериментально. На рис. 7 проведено сравнение депрессии 
упругости паров над системой НГЦ - НЦ, полученное различными авторами, с теорети- 
ческой кривой. 

Из экспериментов известно, что при концентрации НГЦ свыше 0,5 в системе 
НГЦ — НЦ изменение его химического потенциала, по сравнению с чистым вещест- 
вом, практически равно нулю. В то же время теория не дает термодинамического пре- 
дела существования раствора. Возможное объяснение этих факторов следует искать в 
кинетических особенностях образования системы НЦ - растворитель. 

При вводе растворителя в НЦ, когда его концентрация еще достаточно мала, 
практически все его нитратные группы контактируют с нитратными группами НЦ. С 
повышением концентрации растворителя до стехиометрической по функциональным 
группам отношение гетерогенных нитрат-нитратных контактов к числу гомогенных 
стремится к отношению соответствующих термодинамических вероятностей. Дальней- 
шее повышение концентрации растворителя приводит к практически полному насыще- 
нию всех функциональных групп примерно к 45—50% весовой концентрации раство- 
рителя. Следовательно, энергетический фактор, обусловленный взаимодействием гете- 
рогенных контактов, к этому моменту исчерпывается (дифференциальная энтальпия 
растворителя становится равной нулю). Остается лишь энтропийный фактор, который, 



Рис.6. Тепловой эффект растворения 
НГЦ в НЦ 


Рис. 7. Изменение депрессии упруго- 
сти пара НГЦ над системой 
НЦ-НГЦ: 1 - расчет, 2-4 - экс- 
перимент 
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по-видимому, недостаточен, чтобы при температурах изготовления и выдержки, при- 
меняемых в обсуждаемых экспериментах 70 - 80°С, обеспечить дальнейшее про- 
никновение растворителя в жесткий полимер, уже частично им заполненный до полно- 
го выравнивания концентрации. Снятие подогрева, при котором система выдерживает- 
ся, приводит к «замораживанию» незакончившегося процесса проникновения раство- 


рителя. 


О ВОЗМОЖНОСТИ ДИСПЕРСНОГО РАСПРЕДЕЛЕНИЯ РАСТВОРИТЕЛЯ 
В КАПИЛЛЯРАХ И ПОРАХ НЦ 

Выше было указано, что зависимость давления паров растворителя над баллиститным 
порохом могла бы отождествляться с изотермой сорбции капиллярно-пористого тела. 

Если ввести предположение о равновесности системы жидкость - капиллярно- 
пористое тело, а также условие, что заполнение пор жидкостью происходит последова- 
тельно, начиная с пор, обладающих наибольшим капиллярным потенциалом [7] и, в ка- 
честве модели капиллярно-пористого тела выбрать жесткое тело с несквозными непе- 
ресекающимися капиллярами постоянного радиуса, то придем к максимальной величи- 
не контактной энергии, выделяющейся в процессе заполнения пор жидкостью, отне- 
сенной к единице объема, т. е. будет дана верхняя оценка тепловыделению при взаимо- 
действии компонентов. 

Площадь поверхности капилляров, занятых жидкостью, в единице объема 
капиллярно-пористого тела составляет 



где г радиус капилляра, 8 - доля площади поперечного сечения, занятая жидкостью. 

С учетом кривой г(со), полученной на основе -Н- (со), Р однозначно определяет- 
ся как функция от со (здесь со - массовая доля полимера, р 0 , р - давление пара над чис- 
тым растворителем и над капилляром соответственно). 

Можно оценить верхнее значение удельной избыточной энергии контакта рас- 
творителя с поверхностью занятых капилляров, предполагая, что все нитратные группы 
поверхностного слоя НГЦ вступили в контакт с нитратными группами НЦ. На основе 
приведенных выше расчетных данных эта величина составляет - ДЪ = 1 .76 • 1 0' 6 Дж/см 2 

Зависимость выделяемой энергии при капиллярной пропитке НЦ нитроглице- 
рином в функции ю: - ДН = - ДЬ Р показана на рис. 6 штриховой линией. 
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По-видимому, модель капиллярно-пористого тела с жесткими стенками капил- 
ляров не вполне точна. В качестве второго крайнего случая была рассмотрена модель с 
эластичными стенками капилляров в предположении, что реальную систему отобража- 
ет некоторая средняя модель. В этом случае радиусы капилляров растягиваются от ну- 
ля, когда в них нет жидкости - сплошное тело, до г, когда появляется в теле дисперги- 
рованная жидкость. Получено максимальное отклонение - ДН на 5 - 7% в сторону 
уменьшения (рис. 6), по сравнению с первой моделью. 

Очевидно, наблюдаемые в эксперименте тепловые эффекты не могут быть объ- 
яснены выделением энергии за счет поверхностного контакта. Можно, однако, предпо- 
ложить, что часть НГЦ диспергируется до молекул, тогда штриховая кривая на рис. 6 
должна ответвляться от сплошной кривой, отвечающей истинному раствору, при кон- 
центрации растворителя, соответствующей пределу растворимости. Учитывая очень 
пологий характер штриховой линии, приходится сдвинуть точку ответвления практиче- 
ски к точке максимума теоретической кривой, т.е. допустить существование истинного 
раствора, по крайней мере, до концентрации, отвечающей максимуму тепловыделения. 

Таким образом, возможная неоднородность пороховой массы не может быть 
объяснена с точки зрения термодинамической несовместимости компонентов и. по- 
видимому, обязана кинетическим факторам, действующим при изготовлении пороха. 

КИНЕТИКА ПРОНИКНОВЕНИЯ РАСТВОРИТЕЛЯ В НЦ И РАСПРЕДЕЛЕНИЕ 
КОМПОНЕНТОВ В БАЛЛИСТИТНОМ ТОПЛИВЕ 

На основе полученной теоретической зависимости химического потенциала раствори- 
теля от его концентрации в растворе с полимером может ставиться задача о проникно- 
вении растворителя в полимер в процессе изготовления баллиститного топлива. 

Технологический процесс изготовления баллиститных порохов предусматри- 
вает в качестве одного из своих этапов ввод растворителя в размельченную НЦ. пред- 
ставляющую собой рыхлые крупинки размером до 0.2 - 0.5 мм. Можно полагать, что 
скорость набухания такого рыхлоупакованного жесткого полимера, как НЦ, лимитиру- 
ется диффузией растворителя. 

В качестве математической модели используем систему кубов или шаров рых- 
лой структуры, которая изображает размельченную НЦ, с промежутками между ними, 
заполненными растворителем — НГЦ. В начальный момент времени внешняя поверх- 
ность частицы, выделенной в данной системе, вступает в контакт с растворителем ко- 
торый начинает мигрировать в нее, стремясь ликвидировать разность химических по- 
тенциалов. 

В первом приближении положим, что в этом процессе объем частицы не изме- 
няется и количество растворителя отвечает процентному составу топлива. Поэтому до 
некоторого момента времени относительная концентрация растворителя на внешней 
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границе частицы равна 1 , а затем это условие должно быть заменено на условие сим- 
метрии (контакта с соседней частицей). 

Можно представить себе два крайних варианта математической модели систе- 
мы. Первый - в предположении анизотропии частицы НЦ, выраженной в такой степени, 
что растворитель может проникать в нее только с двух противоположных сторон, на- 
пример, для простоты частица в виде куба с ориентацией волокон, перпендикулярной 
двум противоположным граням, через которые только и возможен ввод растворителя. 
Второй вариант - изотропная частица в виде шара с поступлением растворителя по всей 
поверхности сферы. В первом случае задача сводится к задаче о диффузии (теплопро- 
водности) в пластине, вторая - в шаре. 

Уравнение диффузии компонента в неидеальных растворах, когда потенциалом 
переноса служит химический потенциал диффундирующего компонента, имеет сле- 
дующий вид для пластины 


дс _ д ( б Іпа дс л 
ді дх ^ 0 б 1п с дт , 

где с - относительная концентрация компонента; Эо - коэффициенты диффузии в иде- 
альной системе; а - активность диффундирующего компонента; і, г - переменные. 

Активность равна депрессии упругости пара диффундирующего компонента, 
которая берется здесь как функция концентрации по теоретической зависимости давле- 
ния пара растворителя. Аппроксимация этой зависимости дает с достаточной точно- 
стью при всех концентрациях 


бПпа = е к + ьс 
сі\пс 


к = 1.83, Ь = -5.9 


^ С вводом переменных г - г / К, где К. - половина толщины пластины и Ро = 
ЭоУК. -массообменного критерия Фурье, - получаем безразмерную форму уравнения 
(7) и граничные условия 

с(±1 , Ро) — С 0 0< Ро< р0] 

дс 

~ = г (±1,Ро) = 0 РО]<Ро<оо 

дг 


По достижении момента времени Ро, весь растворитель расположенный между 
частицами, полностью расходуется, и каждая частица начинает контактировать, с со- 
седними частицами. Начальное условие 

с(± г , 0) = 0 0 < г < 1 
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На рис. 8 показано распределение концентраций, полученное численно, мето- 
дом конечных разностей. Сплошные линии отвечают Со=1. штрихпунктирные Со=0.5 
Цифрами обозначены значения Ро. 

Расчетному распределению соответствует более глубокое выравнивание кон- 
центраций в центральной области пластины по сравнению с классическим распределе- 
нием (штриховые линии). 

Кроме того, с ростом средней концентрации растворителя в пластине наблю- 
дается более резкое возрастание с но сравнению с классическим случаем, что свиде- 
тельствует о более глубоком замедлении диффузионного процесса. Возможно, указан- 
ным эффектом объясняется поведение термодинамических функций при весовых кон- 
центрациях растворителя свыше 0.5, найденное экспериментально. Скорость растворе- 
ния падает настолько, что полное исчерпание растворителя оказывается практически 
далеким во времени, и система замораживается в неравновесном, двухфазном виде. Та- 
ким образом, условный «предел совместимости» НЦ и растворителя может носить ки- 
нетический характер, если учесть, что термодинамический предел отсутствует. 

Распределение концентраций, полученное при с 0 = 0.5, отражает факт, который 
в классической линейной задаче показался бы парадоксальным - меньшая концентра- 
Ция на границе приводит к более скорому достижению любой величины средней кон- 
центрации и более глубокому выравниванию концентраций по объему по сравнению 
со случаем, где с 0 =1. Это может служить объяснением найденной экспериментально 
закономерности: сорбционный ввод растворителя в НЦ дает большую депрессию упру- 



Рис. 8. Распределение концентрации растворителя в анизотропной частице НЦ 
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гости пара растворителя над топливом по сравнению с капельным вводом, т.е. первый 
способ обеспечивает за то же время более глубокое проникновение и выравнивание 
концентрации растворителя в полимере. 

Аналогично решается задача и для изотропной частицы (задача о диффузии в 
шаре), дающая несколько более глубокое выравнивание концентрации растворителя по 
объему частицы. 

О распределении концентраций по объему баллиститного топлива можно су- 
дить, если привести в соприкосновение все частицы с образованием трехмерной струк- 
туры. Естественно, что при полностью анизотропных частицах в распределение рас- 
творителя по объему вкрадывается элемент случайности, связанный с произвольной 
ориентацией в пространстве каждой частицы. В массе изотропных частиц наблюдается 
упорядоченность в распределении концентраций. Реальное топливо, по-видимому, 
включает в себя частицы разной степени ориентации и разных размеров. 

В определении скорости диффузии растворителя в полимер здесь не учитыва- 
лось изменение реологического противодействия полимера вводу растворителя от кон- 
центрации, учет которого смог бы уточнить полученные результаты. 

ВЫВОДЫ 

Структура исследования, проведенного в 1966 г., состояла в следующем: произведен 
расчет межмолекулярного взаимодействия нитратов и установлен механизм взаимо- 
действия нитратных групп в системе мономер-полимер, в том числе показано, что сте- 
рически невозможно образование максимально прочных связей у НЦ и тепловые эф- 
фекты при растворении НЦ связаны с переходом к более прочным связям НЦ и раство- 
рителя по сравнению с таковыми у чистой НЦ; на базе такого расчета рассмотрено па- 
раллельно две модели НГЦ-НЦ: решеточная модель истинного раствора и модель ка- 
пиллярно-пористого тела, термодинамические характеристики которых служат предме- 
том сравнения с экспериментальными данными, - сделан вывод об отсутствии термо- 
динамических факторов, препятствующих взаимному неограниченному растворению 
компонентов; показано, что предположение о диспергировании растворителя по капил- 
лярам НЦ (двухфазная система) противоречит экспериментальным данным по тепло- 
вым эффектам растворения; рассчитан коэффициент диффузии растворителя в полиме- 
ре как функция относительной концентрации растворителя для модели истинного рас- 
твора, наконец, рассмотрена задача о проникновении растворителя в индивидуальную 
частицу НЦ и сделан вывод о неравномерном распределении компонентов по объему 
топливной массы. Неравномерное распределение компонентов является задающей 
функцией (задающим полем) в реализации локально-нестационарного горения балли- 
ститов. 
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шобе! оі ѵегіІаЫе шопотег-роіушег апб а тосіеі оі саріііагу-ргореііапіз зузіеш аге 
изеё. ТЬе сопсіизюпз із шиіиаі біззоіѵе оГЬазе ргореііапі сошропепіз. ТЬе ргоЫеш оГа 
зоіѵепі зргеаЬ іпіо іпёіѵісіиаі піігосеііиіозе рагіісіе із зоіѵесі. ТЬе сопсіизіоп із ипеѵеп 
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АННОТАЦИЯ. Предложен метод расчета процессов упорядочения и самооргани- 
зации наноструктур с использованием потенциала парного взаимодействия элемен- 
тов наноструктуры, учитывающего силы и моменты сил. Представлены результаты 
исследования самоорганизации системы из двух наночастиц и выполнен анализ ус- 
тойчивости различных типов их равновесия. 


ВВЕДЕНИЕ 

Экспериментальные исследования показывают, что в наноструктурах часто наблюда- 
ются процессы упорядочения и самоорганизации (“зеІР-аззешЫіп^” или “$е1Г- 
ог§апІ 2 Іп§” в англоязычной литературе), приводящие наносистемы к более организо- 
ванной форме [1-5]. В целом, можно выделить три основных типа указанных процес- 
сов: первый - обусловлен формированием регулярной структуры при взаимодействии 
наноструктурных элементов с поверхностью, на которой они расположены; второй - 
является следствием взаимодействия наноструктурных элементов между собой; третий 
- происходит под влиянием среды, окружающей наноструктурные элементы. Влияние 
среды может быть “изотропно распределено’' в пространстве или представлено дейст- 
вием отдельных активных молекул, соединяющих между собой наноэлементы в опре- 
деленном порядке. Кроме того, внешнее воздействие существенно изменяет исходную 
форму структур, образованных наноэлементами. Например, приложение внешнего на- 
пряжения приводит к “вытягиванию” системы наноэлементов по направлению дейст- 
вия максимального растягивающего напряжения; повышение температуры наоборот 
способствует уменьшению пространственной анизотропии наноструктур [2]. Следует 
также отметить, что в процессе самоорганизации наноэлементы, наряду с линейными 
перемещениями, совершают вращательные движения. Последнее можно объяснить 
действием моментов сил, обусловленных несимметрией силовых полей взаимодействия 
наноэлементов, наличием на поверхности наноэлемента локальных областей 
притяжения или “отталкивания”, а также “неизотропным” действием окружающей 
среды. 
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Указанные явления играют важную роль в нанотехнологических процессах. Они 
позволяют разрабатывать нанотехнологии формирования наноструктур методом само- 
сборки (в основе которой находятся процессы самоорганизации) и создавать сложные 
пространственные наноструктуры, состоящие из различных наноэлементов (наноча- 
стиц, нанотрубок, фуллеренов, супермолекул и т.д.) [6]. С другой стороны, стремление 
к самоорганизации, в ряде случаев, препятствует формированию нужной нанострукту- 
ры. Так, наноструктура, возникшая в процессе самоорганизации, как правило, является 
“жесткой” и устойчивой к внешним воздействиям. Например, “слипание” наночастиц 
препятствует использованию изолированных наночастиц в различных нанотехнологи- 
ческих процессах, равномерному перемешиванию наночастиц из разных материалов и 
формированию нанокомпозита с требуемыми свойствами. В этой связи высока акту- 
альность моделирования процессов статического и динамического взаимодействия 
элементов наноструктур. Важным при этом является учет моментов сил взаимодейст- 
вия элементов наноструктур, вызывающих взаимный поворот наноэлементов. 

Задачами данной работы, в развитие ранее выполненных работ автора [7-11], яв- 
лялись постановка задачи расчета процессов упорядочения и самоорганизации наност- 
руктур и анализ их устойчивости. В работе использовано понятие потенциала взаимо- 
действия элементов наноструктуры, дополненное учетом моментов сил взаимодействия 
наноэлементов. 

ПОСТАНОВКА ЗАДАЧИ 

Задача упорядочения и самоорганизации элементов наноструктуры включает три ос- 
новных этапа: первый - построение внутренней структуры и равновесной конфигура- 
ции (формы) каждого отдельного, свободного от воздействий элемента наноструктуры; 
второй - расчет парного взаимодействия двух элементов наноструктуры; третий - уста- 
новление закономерностей пространственного строения и эволюции во времени нано- 
структуры в целом. Рассмотрим подробно формулировку задачи для каждого из ука- 
занных этапов. 

Расчет внутренней структуры и равновесной конфигурации наноэлементов 

Решение задачи расчета структуры и равновесной конфигурации наноэлемента выпол- 
ним методом молекулярной динамики, рассматривая взаимодействие всех образующих 
элемент атомов. Так как на первом этапе решения на наноэлемент не действуют внеш- 
ние силы, то с течением времени он принимают равновесную конфигурацию, которая 
далее используется для следующего этапа расчетов. 

На данном этапе движение атомов, входящих в наноэлемент, определяется систе- 
мой дифференциальных уравнений Ланжевена [12]: 
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N к 

пи- Х і = Т.Г.і+Г,(‘)-а,Щ ' , >' = 1,2,.., ЛГ 4 , 

7=1 ™ 


СІІ 


( 1 ) 


при краевых условиях 


х , = х ю , X - Ѵ і0 . 1 = 0, х, с= О к , (2) 

где Ы к - число атомов, составляющих каждую наночастицу; ггц- масса і-го атома; 
Х І 0 ,Х; - начальные и текущие координаты і-го атома, соответственно; К - силы меж- 
атомного взаимодействия; О к - область, занимаемая соответствующей наночастицей; 
Ѵ і0 , V - начальная и текущая скорости і-го атома, соответственно; ОС - коэффициент 

'‘трения” в атомарной системе; Р (і)- случайный набор сил при расчетной температуре, 
задаваемый распределением Гаусса. 

Силы межатомного взаимодействия, обычно, являются потенциальными и опре- 
деляются из соотношения 


^ ? фГ •' Х "- і 1,2 к ' 


( 3 ) 


где Ру - радиус-вектор, определяющий положение і-го атома относительно і-го атома; 

Ф(Ру ) ‘ потенциал, зависящий от взаимного расположения всех атомов; п - число ти- 
пов взаимодействия между атомами. 

Потенциал Ф( ру ) , в общем случае, задается в виде суммы нескольких состав- 
ляющих, соответствующих различным типам взаимодействия: 

Ф(Ріі) = ФсЬ + Ф., + Ф,а + Фр 8 + Ф„ + Ф еі + Ф„ь- (4) 

Здесь имеются в виду потенциалы: Ф сЬ - химических связей; Ф ѵа - валентных углов; 
Фіа" торсионных углов; Ф рё - плоских групп; Ф ѵѵ - Ван-дер-ваальсовых контактов; 

Ф е5 - электростатики; Ф ьь - водородных связей. 

Указанные слагаемые имеют различный функциональный вид, а значения пара- 
метров для потенциалов взаимодействия определяются на основе экспериментов (кри- 
сталлографических, спектральных идр.) и квантомеханических расчетов [12-13]. 
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Расчет парного взаимодействия двух элементов наноструктуры 

На данном этапе решения задачи рассматриваем два взаимодействующих наноэлемен- 
та, расположенные друг от друга на расстоянии 8 С (рис.1) между их центрами масс и 
ориентированных относительно друг друга под определенными заданными углами. 

Взаимодействие атомов, входящих в наноэлементы, приводит к относительному 
перемещению и повороту последних. Соответственно, сумма всех сил взаимодействия 

атомов наноэлементов, в общем случае, приводится к главному вектору сил Р и 

главному моменту М . 

- - ~ ". " 2 

Г =Гы=-Рь 2=^ (7) 

'=1 7=1 


X 




Рис.1. Два взаимодействующих наноэлемента; М, Р - главный момент и главный 
вектор сил, соответственно 
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М = Мь\ = -Мы , 


( 8 ) 


где і, ] - атомы первой и второй наночастицы, соответственно. 

Главной целью данного этапа расчета является построение зависимостей сил и 
моментов взаимодействия элементов наноструктуры от расстояния между центрами 
масс элементов наноструктуры 8 С , углов взаимной ориентации наноэлементов 

Ѳ,,Ѳ 2 ,Ѳ 3 (формы наноэлемента) и характерного размера наноэлемента I) . В общем 
случае эти зависимости можно представить в виде 



( 8 ) 



( 9 ) 


Для сферических наноэлементов углы взаимной ориентации не влияют на силу их 
взаимодействия, поэтому в уравнении (9) момент тождественно равен нулю. 

В частности, расчеты данной зависимости для сферических наночастиц, приве- 
денные в [11] показали, что в этом случае функция Р ы может быть записана как про- 
изведение функций, зависящих от размера наночастиц и расстояния между ними: 


^=ф,і№)-ф п ( О ) 


(10) 


В общем случае функции (8) и (9) могут быть аппроксимированы, по аналогии с 
(10), как произведение функций ^.©„©^©з,/?, соответственно. Следует отметить, 
что, учитывая дальнейшее численное решение задачи самоорганизации наноэлементов, 
достаточно задать указанные функции в табличном виде и использовать линейную (или 
нелинейную) интерполяцию их в пространстве. 

Расчет процессов упорядочения и самоорганизации наноструктуры 

При исследовании процессов упорядочения и самоорганизации элементов нанострук- 
туры рассматриваем движение каждого элемента системы, как единого целого. При 
этом поступательной движение центра масс каждого наноэлемента задаем в системе 
координат Хі, Х 2 ,Хз, а вращение в наноэлемента описываем в системе координат 1\, Ъ 2 
2з, связанной с центром масс наноэлемента (рис.1). Система уравнений, описывающая 
данные процессы, имеет вид 
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где Х і ,Ѳ; - координаты центров масс и углы ориентации в пространстве главных осей 
2 Ь 2 2 ,2 3 наноэлементов; - силы взаимодействия наноэлементов; 

- внешние силы, действующие на наноэлементы; N. - число наноэле- 
ментов, М ^ - масса наноэлемента; А4 Л/ ^ * момент сил взаимодействия нано- 

элементов; М 2 ^М ' 2 - внешние, воздействующие на наноэлементы, моменты; 
^ 2^>^ 2^>^ 2 ~ моментз инерции наноэлемента. 

Начальные условия для системы уравнений (1 1), (12) имеют вид 


= ®/о> Ѵ( = Ѵ ю , Ѳ, = Ѳ і0 , і = 0. 


( 13 ) 
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ПРИМЕРЫ И АНАЛИЗ РЕЗУЛЬТАТОВ РАСЧЕТОВ 

Реализацию предложенной выше методики рассмотрим на примере расчета самоорга- 
низации системы из двух кубических наночастиц, взаимодействие атомов которых оп- 
ределяется потенциалом Морзе[12]. На рис. 2 показаны возможные взаимные положе- 
ния этих наночастиц. Позиции, в которых главный момент сил равен нулю, соответст- 
вуют парам наночастиц 2-3; 3-4; 2-5 (рис.2) и определяют возможные положения их 
равновесия. 

На рис.З приведена зависимость момента силы взаимодействия между кубиче- 
скими частицами 1-3 (рис. 2) от угла их относительного поворота. Из графика следует, 
что при угле поворота частицы 1 относительно частицы 3, равном л/4, момент сил их 
взаимодействия равен нулю. При увеличении или уменьшении данного угла появляется 
момент сил. В диапазоне л/8<Ѳ<Зл/4 этот момент мал. Вне данного диапазона момент 
сил быстро возрастает. Следует особо отметить существенную роль в установлении 
равновесия наночастиц расстояния 8 между ними. Если 8>8 0 (где 8 0 — расстояние, на 
котором силы взаимодействия между наночастицами равны нулю), частицы притяги- 
ваются друг к другу. При этом знак момента соответствует знаку отклонения угла Ѳ от 
л/4. При 8<8 0 (отталкивание наночастиц) знак момента противоположен знаку откло- 
нения угла. 1 о есть, в первом случае увеличение отклонения угла вызывает увеличение 
момента, обусловливающего движение наноэлемента в данном направлении, а во вто- 
ром случае отклонение угла вызывает увеличение момента, препятствующего движе- 
нию наноэлемента в данном направлении. Таким образом, первый случай соответствует 



Рис. 2. Характерные позиции кубических наночастиц 
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Рис. 3. Зависимость момента силы М(нм) взаимодействия между кубическими 
частицами 1-3 (показаны на рис. 2) от угла их относительного поворота 
Ѳ(рад) 

неустойчивому равновесию, а второй устойчивому равновесию наночастиц. Влияние 
параметра 3 при этом иллюстрируют графики изменения потенциальной энергии сис- 
темы взаимодействия двух кубических наночастиц (рис.4). Здесь кривая 1 соответству- 
ет условию 5<8 0 и имеет хорошо выраженный минимум в области 0.3<Ѳ < 1.3. При Ѳ 
< 0.3 и Ѳ > 1.3 потенциальная энергия взаимодействия резко возрастает, что приводит к 
возвращению системы в исходное равновесное положение. При 8>8 0 ( кривые 2-5) 

график потенциальной энергии в точке Ѳ = 0 имеет максимум, что соответствует е не- 
устойчивому положению. 

Проведенный теоретический анализ подтверждается работами ученых из универ- 
ситета штата Нью Джерси и университета Калифорнии в Беркли, которые эксперимен- 
тально обнаружили самоорганизацию кубических микрочастиц цирконата - титаната 
свинца (Р2Т) [14]: упорядоченные коллективы кубических микрокристаллов из гидро- 
термально-синтезированного Р2Т сформировали на границе раздела воздух-вода пло- 
ский слой частиц, причем частицы заняли наиболее устойчивое положение, соответст- 
вующее позиции 2-3 на рис.2. 

Таким образом, анализ взаимодействия двух кубических наночастиц показал, что 
возможны различные варианты их конечного стационарного состояния равновесия, в 
котором главные векторы сил и моментов равны нулю. Однако, существуют как устой- 
чивые, так и неустойчивые стационарные состояния данной системы: положение час- 


226 


ХИМИЧЕСКАЯ ФИЗИКА И МЕЗОСКОПИЯ. Том 7, № 2 


МОДЕЛИРОВАНИЕ ПРОЦЕССОВ УПОРЯДОЧЕНИЯ И САМООРГАНИЗАЦИИ НАНОСТРУКТУР 



Рис.4. Графики изменения потенциальной энергии Е (нм) взаимодействия двух 
кубических наночастиц от угла их относительного поворота Ѳ (рад) и рас- 
стояния между ними (положение наночастиц 1-3, рис. 2) 

тиц 2-3 -устойчивое, а позиции 3-4, 2-5 обладают ограниченной устойчивостью или не- 
устойчивы, в зависимости от расстояния между наночастицами. 

Следует отметить, что для структур, состоящих из большого количества наноча- 
стиц, может существовать множество устойчивых стационарных и неустойчивых форм 
равновесия. Соответственно, могут возникать устойчивые и неустойчивые нанострук- 
туры композиционных материалов. Поиск и анализ параметров, определяющих форми- 
рование устойчивых наносистем, является актуальной задачей. 
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РЕФЕРАТЫ 


УДК 623.459.454.0004.74 

Влияние степени извлечения сортируемых фракций на теплотехнические свойст- 
ва бытового мусора. Петров В.Г., Чечина А. А. ХИМИЧЕСКАЯ ФИЗИКА И 
МЕЗОСКОПИЯ. 2005. Том 7, № 2. С. 135-145. 


На примере г. Ижевска показано влияние извлечения сортируемых фракций на 
теплотехнические характеристики твердых бытовых отходов (ТБО). Оценено содержа- 
ние в ТБО влаги, горючих веществ и золы в разные сезоны года, а также их морфоло- 
гический состав. Приведены результаты расчетов теплоты сгорания ТБО. В результате 
сортировки происходит изменение свойств получаемого субстрата ТБО, что может ска- 
заться при их утилизации различными способами. Табл. 3. Ил. 2. Библиогр. 9. 

УДК 532.613.5 

Метод расчета коэффициента поверхностного натяжения. Корепанов М.А. 
ХИМИЧЕСКАЯ ФИЗИКА И МЕЗОСКОПИЯ. 2005. Том 7, № 2. С. 146-154. 

Предлагается метод расчета коэффициента поверхностного натяжения через те- 
плофизические свойства вещества. Предлагаемый метод позволяет рассчитывать коэф- 
фициент поверхностного натяжения как неполярных, так и полярных жидкостей через 
У чет дипольного момента их молекул. Показано хорошее совпадение результатов рас- 
четов коэффициента поверхностного натяжения жидкостей с экспериментальными 
данными. Табл. 3. Библиогр. 1 1. 

УДК 544.344.01 

Расчет термодинамических свойств метастабильных конденсированных класте- 
ров МЬ 2 -1ѴЪ5 и состава жидкого ниобия с учетом их существования Моисеев Г.К. 
ХИМИЧЕСКАЯ ФИЗИКА И МЕЗОСКОПИЯ. 2005. Том 7, № 2. С. 155-161. 

С применением данных молекулярной динамики и известных расчетных мето- 
дов оценены термодинамические свойства метастабильных конденсированных класте- 
ров N62 - ХЬ 5 ; методами термодинамического моделирования с учетом атомов и кла- 
стеров рассчитан состав жидкого ниобия при Р = 9.8066- ІО' 2 МПа иТ = 2800-5300 К: Т 
и АН кипения. Ил. 1. Табл. 3. Библиогр. 13. 
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УДК 621.385.833 

Исследование области разрыва «шейки» заготовки зондирующей иглы СТМ при 
ее изготовлении методом химического травления. Липанов А.М., Шелковников 
Е Ю., Тюриков А. В., Д.В.Гудцов. ХИМИЧЕСКАЯ ФИЗИКА И МЕЗОСКОПИЯ. 2005. 
Том 7, №2. С. 162-168. 

В статье рассмотрена полученная с применением метода молекулярной динами- 
ки модель формирования острия зондирующей иглы СТМ. Приведены результаты ис- 
следования влияния силы тяжести нижней части заготовки иглы на характер разрыва ее 
«шейки» с целью получения атомарных микровыступов в области разрыва. Ил. 3. Биб- 
лиогр. 9. 

УДК 622.276.72 

Магнетизм нефти и омагничивание ее для уменьшения асфальтосмолопарафино- 
вых отложений. Бородин В.И., Хрущев А.Д., Зинин А. В., Тарасов Е.Н. 
ХИМИЧЕСКАЯ ФИЗИКА И МЕЗОСКОПИЯ. 2005. Том 7, № 2. С. 169-180. 

Приведены результаты исследований воздействия постоянного магнитного поля 
на флюид нефть плюс вода с целью построения эффективных магнитных устройств. 
Показано, что применение магнитной обработки позволяет существенно уменьшить 
образование отложений на стенках насосно-компрессорных труб, и увеличить срок ра- 
боты скважин. Ил. 7. Библиогр. 12. 

УДК 662.615.1 

О пространственной компонентной неоднородности базовых баллиститных ком- 
позиций в аспекте ее влияния на скорость горения. Романов О.Я. ХИМИЧЕСКАЯ 
ФИЗИКА И МЕЗОСКОПИЯ. 2005. Том 7, № 2. С. 201-218. 


Осуществлялась проверка гипотезы влияния неравномерного распределения по 
объему основных компонентов баллиститов, которые обычно считаются гомогенными 
композициями, на реализацию локально-нестационарного механизма горения. Произ- 
веден расчет межмолекулярного взаимодействия нитратов, рассмотрена решеточная 
модель истинного раствора мономер-полимер и модель капиллярно-пористого тела. 

Сделан вывод об отсутствии термодинамических факторов, препятствующих 
взаимному неограниченному растворению компонентов. Решена задача о проникнове- 
нии растворителя в индивидуальную частицу НЦ. Сделан вывод: неравномерное рас- 
пределение компонентов является задающей функцией (задающим полем) в реализации 
локально-нестационарного горения баллиститного топлива. Табл. 1. Ил. 8. Библиогр. 7. 
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УДК 539.3: 519.68 


Моделирование процессов упорядочения и самоорганизации наноструктур. 

Вахрушев А.В. ХИМИЧЕСКАЯ ФИЗИКА И МЕЗОСКОПИЯ. 2005. Т.7. №2. С.219-228. 

Предложен метод расчета процессов упорядочения и самоорганизации наност- 
руктур с использованием потенциала парного взаимодействия элементов нанострукту- 
ры, учитывающего силы и моменты сил. Представлены результаты исследования само- 
организации системы из двух наночастиц и выполнен анализ устойчивости различных 
типов их равновесия. Ил. 4. Библиогр. 14. 
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